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We present a general framework for understanding and an- 
alyzing critical behaviour in gravitational collapse. We adopt 
the method of renormalization group, which has the following 
advantages. (1) It provides a natural explanation for various 
types of universality and scaling observed in numerical stud- 
ies. In particular, universality in initial data space and univer- 
sality for different models are understood in a unified way. (2) 
It enables us to perform a detailed analysis of time evolution 
beyond linear perturbation, by providing rigorous controls on 
nonlinear terms. Under physically reasonable assumptions we 
prove: (1) Uniqueness of the relevant mode around a fixed 
point implies universality in initial data space. (2) The criti- 
cal exponent /3bh and the unique positive eigenvalue k of the 
relevant mode is exactly related by /3bh = /3/k, where /3 is a 
scaling exponent used in calculating the eigenvalue. (3) The 
above (1) and (2) hold also for discretely self-similar case (re- 
placing "fixed point" with "limit cycle"). (4) Universality for 
different models holds under a certain condition. [These are 
summarized as Theorems |l[ |^ and ^ of Sec. |^.] 

According to the framework, we carry out a rather com- 
plete (though not mathematically rigorous) analysis for per- 
fect fluids with pressure proportional to density, in a wide 
range of the adiabatic index 7. The uniqueness of the rele- 
vant mode around a fixed point is established by Lyapunov 
analyses. This shows that the critical phenomena occurs not 
only for the radiation fiuid but also for perfect fiuids with 
1 < 7 < 1.88. The accurate values of critical exponents are 
calculated for the models. In particular, the exponent for the 
radiation fluid /3bh — 0.35580192 is also in agreement with 
that obtained in numerical simulation. 

PACS numbers: 04.40.-b, 04.70.Bw, 64.60.Ak 



I. INTRODUCTION 

Gravitational collapse with formation of black holes is 
one of the main problems of classical general relativity. 
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A very important question here is how formation of a 
black hole depends on initial data, for example, which 
initial data evolve into black holes and which do not. 
It is concerned with many important problems such as 
cosmic censorship [|l||. 

Christodoulou Pppl rigorously analyzed such a prob- 
lem in a relatively simple system, namely, spherically 
symmetric scalar field collapse. He proved that the field 
disperses at later times if the initial field is sufficiently 
weak, and also found a sufficient condition for black hole 
formation. However, the behaviour of the field and the 
space-time for initial data sets around the threshold of 
the black hole formation remained open, and he posed 
a question whether the phase transition between space- 
times with and without black holes is of the first or of 
the second order, i.e., whether the mass of the black hole 
formed is a continuous function at the threshold. 

Choptuik dl numerically analyzed the system and ob- 
tained a result which strongly suggests the latter. Fur- 
thermore, he discovered a "critical behaviour" of the 
gravitational collapse of the system. His result can be 
summarized as follows: Let the initial distribution of the 
scalar field be parametrized smoothly by a parameter p, 
such that the solutions with the initial data p > Pc con- 
tain a black hole while those with p < pc do not. For 
several one-parameter families investigated, near-critical 
solutions (p Ri Pc) satisfy the following: (1) the critical 
solution (i.e. p = Pc) is universal in the sense that it 
approaches the identical space-time for all families, (2) 
the critical solution has a discrete self-similarity, and (3) 
for supercritical solutions {p > Pc) the black hole mass 
satisfies Mbh oc (p — Pc)^^^ and the critical exponent 
Pbh, which is about 0.37, is universal for all families. 
Abrahams and Evans Q found similar phenomena in ax- 
isymmetric collapse of gravitational wave and obtained 
/?BH = 0.38. Evans and Coleman Q found similar phe- 
nomena with /3bh — 0.36 in spherically symmetric col- 
lapse of radiation fluid, in which case the self-similarity 
is not discrete but continuous. Employing a self-similar 
ansatz they also found a numerical solution which fits the 
inner region of the near-critical solutions very well. Since 
the above values of /3bh were close to each other, some 
people considered that there is a universality over many 
systems. 
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For a deeper understanding of these phenomena, an 
important direction is to pursue the analogy not only 
phenomenologically but also theoretically between the 
"critical phenomena" in gravitational collapse and those 
in phase transitions in statistical mechanics. Argyres |^] 
reviewed the above numerical results in terms of rcnor- 
malization group (RG). Independently, the present au- 
thors ||ic| ] showed that both qualitative and quantitative 
understanding is possible by introducing renormalization 
group ideas. They showed that uniqueness of the relevant 
(unstable) eigenmode of the linear perturbations around 
the critical self-similar space-time is essential for the uni- 
versality to be observed and that the critical exponent 
Pbh is given by /3bh = where k is the only eigenvalue 
with positive real part. Emphasized there was that the 
standard renormalization group argument implies that 
the linear perturbation analysis is sufficient to obtain the 
exact value of /3bh- They confirmed that their scenario 
holds in the case of radiation fluid and obtained a very 
accurate value of the exponent, /3bh — 0.35580192. The 
difference from those obtained for scalar field collapse 
(0.37) and gravitational wave collapse (0.38) was beyond 
the possible numerical errors in their simulations. This 
showed that there is no universality between radiation 
fluid collapse and other systems. Maison applied 
the method to the systems of perfect fluids with pressure 
proportional to density and showed how /3bh depends 
on the equation of state, under the assumption that the 
same phenomena occur and the same scenario holds as 
the radiation fluid. Gundlach |l2j confirmed /3bh = ^/k 
in scalar field collapse. 

Some new and intriguing phenomena in the critical 
behaviour have been found such as bifurcation [ p^pl) 
and coexistence of first and second order phase tran- 
sitions [ p^ . However, in this paper we will pursue a 
slightly different direction, namely concentrate on pro- 
viding deeper understanding of basic aspects of the crit- 
ical behaviour. 

This paper is the fully expanded version of Ref. [p^ . 
(Some of the materials here have already been treated in 
|6|.) The aim of the paper is twofold. First, we present 
a general, mathematically rigorous, framework in which 
critical behaviour observed in various models so far can 
be analyzed and understood in a unified way, based on 
a renormalization group picture. Now that there seems 
to be some consensus on the qualitative picture on the 
simplest case, we in this paper focus our attention to 
two problems which remain to be more clarified. (1) We 
give a sufficient condition under which the critical expo- 
nent is given exactly by /3bh — P/k. In the above, /3 
is a scaling exponent which appears in calculating the 
Lyapunov exponent (see Sec. |l for details). Note that in 
pO| we only considered the radiation fluid collapse, where 
one can take (3 = 1. (2) We also consider the problem 
of universality, in the wider sense of the word that de- 
tails of the definition of the model is irrelevant for the 
critical behaviour, and we present a sufficient condition 
under which the universality of this kind occurs. Both 



conditions are sufficiently weak that they are expected 
to be applicable to most systems for which the critical 
behaviour has been observed. 

Second, we present a more thorough analysis of the 
perfect fluid collapse, with various adiabatic index 7. In 
addition to presenting the details of our analysis on the 
radiation fluid collapse reported in [ p^ , we here report 
a result of a new method of analysis, the so called Lya- 
punov analysis, which was performed to further conflrm 
the uniqueness of the relevant mode around the critical 
solution. A Lyapunov analysis and the shooting method 
for an ordinary differential equation adopted in ]T^ ] are 
complementary to each other in the following sense: (1) 
The former extracts eigenmodes in the descending order 
of its real part, whereas the latter affords information 
only of a finite region of complex n plane. (2) The Lya- 
punov analysis is useful even if the eigenmodes do not 
form a complete set. (3) It is easier in the latter than 
in the former to numerically obtain accurate values of 
the eigenvalues n hence the critical exponent /3bh- Our 
analysis in the second part is sufficiently convincing for 
physicists, we believe. But mathematically it is unsat- 
isfactory in the sense that it does not provide rigorous 
theorems. Rigorous full analysis of the system is a sub- 
ject of future studies. 

This paper is organized as follows. First, we present 
our general scenario in Sec. We define renormaliza- 
tion group transformations acting on the phase space 
(Sec. HA), and see how all the aspects of critical phenom- 
ena for radiation fluid collapse, etc. are deduced from the 
behaviour of the RG flow near a fixed point (Sec. [II B| ). 
We present a sufficient con dition for the exactness of the 
relation /3bh = P/n in Sec. [IC. We then consider more 
general cases, i.e. discrete self-similarity (Sec. [ID) and 
problem of universality for different systems (Sec. HE). 

Then we turn our attention to detailed study of perfect 
fluid collapse. After reviewing the equations of motion 
in Sec. Ill, we first study self-similar solution in Sec. |rv|. 



Then we study perturbations around the self-similar so- 
lution and confirm our picture to an extent in Sec. ^ 
by numerical study. In Sec. VI we establish uniqueness 
of the relevant mode by the Lyapunov analysis. The 
uniqueness of the relevant mode implies that that self- 
similar solution i s res ponsible for the critical behaviour. 
So we list in Sec. VI] the value of the critical exponent /3 



for various 7. Sec. VIII is for conclusions and discussions. 



II. SCENARIO BASED ON RENORMALIZATION 
GROUP IDEAS 

In this section we present a scenario of critical be- 
haviour based on renormalization group ideas, which 
gives clear understanding of every aspect of the phenom- 
ena. We first introduce renormalization group transfor- 
mations acting on the phase space of partial differential 
equations (PDEs) in Sec. II A . In Sec. p^IB| we give a (non- 
rigorous) scenario for the simplest case where the system 



2 



has scale invariance and a (continuously) self-similar so- 
lution plays an important role. Even in this simplest 
case we can argue most of the essential features of the 
critical behaviour. This case covers a perfect fluid col- 
lapse treated in det ail in the subsequent sections. We 
then proceed in Sec. II C to make our scenario into a rig- 
orous mathematical framework which starts from basic 
assumptions, based on renormalization group ideas. The 
case where discrete self-similar solution becomes relevant 
such as scalar field collapse is underst ood e ssentially in 
the same way, as is explained in Sec. [ID, In general 



cases where the system is not scale invariant, the renor- 
malization group transformation drives the equations of 
motion to a fixed point. Renormalization group method 
then tells us which systems can exhibit t he s ame critical 
behaviour, as will be explained in Sec. HE. As exam- 



ples, we consider perfect fluid with a modified equation 
of state, and also show that scalar fields with any poten- 
tial term (which depends on the scalar field only) should 
show the same critical behaviour. 

Currently, there seem to be two kinds of renormaliza- 
tion group type analysis of PDEs, |l7| and [Q. Although 
these two are certainly philosophically related, they dif- 
fer in details. We here adapt the approach of [Q, and 
extend their work to cases where the critical behaviour 
is observed. It would also be very interesting to consider 
applications of the method of n% . 



A. Renormalization group transformation 

We give a general formalism which deals with time evo- 
lution of initial data considered as a flow of a renormal- 
ization group transformation. The argument is general 
but the notations are so chosen in order to suit the case 
considered in the subsequent sections. In particular, we 
explicitly consider only the case of spherical symmetry, 
although our framework can in principle be applied to 
more general cases in n-space dimensions with suitable 
changes: The independent space variable r is increased 
to ri , r2 , . . . , r„ , and the differential operators TZ and T 
become partial differential operators on a space of func- 
tions of n variables. 

We are interested in the time evolution of n unknown 
functions u = (ui, U2, u„), which are real-valued func- 
tions of time t and spatial coordinate r, and which satisfy 
a system of partial differential equations (PDEs) : 



, du du 



0. 



(2.1) 



This real formulation does not restrict applicable sys- 
tems in any way, because any PDEs for complex unknown 
functions can be rewritten in the real form by doubling 
the number of components. We employ this real form, 
because this makes the dimensional counting (e.g. num- 
ber of relevant modes) si mpler a nd clearer. See a Re- 
mark after Claim 1 of Sec. II B 1. (However, in practical 



problems, e.g. complex scalar fields, one could proceed 
without employing the real formulation, as long as some 
caution is taken in counting dimensions.) 

The time evolution of this system is determined once 
one specifies the values of u{t, •) at the initial time t = —1. 
We call the space of functions u(i, •) (at fixed time t) the 
phase space F, which is considered as a real vector space. 
We leave the definition of F (e.g. how smooth functions in 
F should be) rather vague, until we d evelop more detailed 



theory of time evolution in Sec. II B 



To introduce a renormalization group transformation, 
we first define a scaling transformation S{s,a,(3) which 
depends on real parameters s,a = (ai, a2, ■ ■ ■ , Q;„), and 

5(s,a,/3) : u,{t,r) ^ u["^"-'^\t,r) = e"'"u,(e-"i, e'/^V) 

(2.2) 

(As will be clear /3 here is related with, but not equal to, 
/3bh.) 

Now we make the fundamental 



Assumption S ( Invariance under scaling) The sys- 
tem of PD Es (|2.l| ) are invariant under the scaling trans- 
formation ( p. 21) in the following sense: If w is a solution of 
(2.1), then for suitably chosen constants a, (3, and for an 
arbitrary scale parameter s (close to 0), the new function 



u 



(s,Q,/3) 



is also a solution of the system (UA 



In the following we denote by a, /? the constants which 
mak e the above assumption hold for the system of PDEs 
( ^.l|) , and denote u^*'"''') simply by u'*-* unless otherwise 
stated. We consider only those systems for which the 
above Assu mption S holds in Sec. II A through Sec. [ID. 
In Sec. HE, we consider systems for which Assumption S 
does not hold, and address the problem of universality in 
its wider sense. 

We now define a renormalization group transformation 
(RGT) acting on the phase space F as 



K.t/f)(e)=u|^)(-l,e) = e"^^«,(-e- 



-^■'0, (2.3) 



where a, (3 are constants of Assumption S, s > 0, and u 
is a solution of the PDE system. That is, is given by 
developing the initial data u(— 1, r) = U{r) from t = —1 
to t = — e"** by the PDE, and rescaling the spatial coor- 
dinate r and the unknown functions Ui. We have written 
the argument of C/ as ^ rather than r, in order to draw 
attention to the fact that the physical (or geometrical) 
length is given by 



(2.4) 



rather than by ^ itself. 

Because of our assumption of the invariance of the 
PDEs under scaling, TZs forms a (semi) group: 
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7^ 



S1+S2 



Tls2 ° T^Si- 



(2.5) 



From this it follows immediately that the s-derivative of 
TZs at s = is in fact an infinitesimal generator of TZs : 



n = lim 



7^, - 1 



TZs = cxp 



(2.6) 



A fixed point of the renormalization group {TZs \ s G 
M} is a point U* in T satisfying UsiJJ*) = U* for all 
s > 0, and can also be characterized by 'R.{U*) = 0. A 
function u is called self-similar with parameters (a, /3) if 



it satisfies u{t, r) 



{t,r). Each self-similar so- 



lution Uss of the PDEs with parameters (a, P) of As- 
sumption S is related to a fixed point U* of TZs by 
Ussit,r) = {-t)"U*{r{~t)-^). 

The tangent map (or Frechet derivative) Ts^jj of TZs at 
[/ is defined by 



Ts,u{F) = lim 

£-►0 



TZs{U + eF) - TZsU 



Note that T^.t/ is a linear operator for fixed s and ?7, 
in contrast with TZs- The above definition suggests the 
following formal relation: 



TZs{U + F) =TZs{U) +Tsm{F) + 0{Ff 



(2.8) 



where 0{F)'^ denotes a term whose norm is of order 
(with some suitable norm || • ||). 
We define the s-derivative ofTg^u at U as 



Ts 



s.U 



1 



71/ = lim 

which is formally related with TZ by 

TZ{U + F)^ TZ{U) + fu{F) + 0{Ff. 



At a fixed point U* of TZs, ( |2.7] ) simplifies to 
T.(F).lim^-(^^+-^)-^\ 



(2.9) 



(2.10) 



(2.11) 



and we denote the tangent map at U* simply by Ts in the 
following. Ts (at U*) forms a (semi)group like TZs, and 
can be represented in terms of its infinitesimal generator 
f as 



T = Ttj, 



lim 



T^l 



Ts = exp {sT^ . (2.12) 



An eigenmode F oiT with eigenvalue k is a (complex- 
valued) function satisfying (k G C) 



TF ^ kF. 



(2.13) 



Note that we here consider eigenmode (and spectral) 
problem on complexification of F, i.e. we allow complex 
eigenvalues and eigenvectors. This is necessary for our 



later analysis. Although complex eigenvectors are not el- 
ements of F, we can express vectors in F by taking the 
real part of a linear combination of complex eigenvectors. 

These modes determine the flow of the RGT near the 
fixed point U*^. A mode with Re k > 0, a relevant mode, 
is tangent to a flow diverging from U*, and one with 
Re K < 0, an irrelevant mode, to a flow converging to it. 
A Re K — mode is called marginal. 

Remark. Many equations of gravitational systems, in- 
cluding those of the perfect fluid, are not exactly of the 
form of (|]^). More precisely, we have two kinds of un- 
knowns u = {ui,U2, ■ ■ ■ ,Un) and c — (ci, C2, . . . , c™), 
which satisfy the following PDEs of two different classes: 
(l)the evolution equations: 



du du dc . 
ot or or 



(2.7) and (2) the constraint equations: 



Lc { u,c,^,r ] = 0. 



(2.14) 



(2.15) 



In this case, the well posed Cauchy problem is to spec- 
ify the initial values of Ui{—1, •) for i = 1,2, ... ,n (and 
appropriate boundary conditions for q's). In this sense, 
the phase space variables are Mi's, while c^'s should be 
considered as constraint variables. 



Solving ( 2.15 ) for constraint variables Cj's in te rms of 
integrals of u^'s and substituting them intoJ 2.14 ) yields 
evolution equations for m^'s of the type of ( |2.1| ), with the 
difference that u in the argument of L now represents 
some complicated functionals (integrals) of u^'s. In this 
form, we can apply our formalism of the renormalization 
group transformation. 

Remark. In the above, we have presented the formalism 
in (t, r) coordinate with unknowns Ui, in order to directly 
reflect our physical interpretation. However, mathemat- 
ical treatment can in principle be made much simpler by 
introducing new variables 

X = Inr — p\Ti(—t), s = — ln(— t), 

v,{s,x) = e°''''u,{t,r). (2.16) 

In terms of these, the scaling transformation S{s',a,(3) 
becomes a simple translation in s: 



Vi{s,x) 1-^ \s,x) = Vi{s + s',x), 



(2.17) 



^Mathematically, in order to determine the behaviour of the 
flow around U* , we have to know not only the point (dis- 
crete) spectrum but al so th e continuous spectrum of T. See 
Assumption LI of Sec. [IC for details. 



4 



and the parameters a, (3 disappear. One can thus identify 
the scale parameter s' with the new "time" coordinate s. 
The assumed scale invariance of the system implies that 
the EOM do not include s explicitly. Now TZ and T can 
be expressed by usual s-derivatives on the functions Vi . In 
Ip^ , we formulated our rcnormalization group approach 
directly in terms of these variables. 



B. Scenario 

In this subsection we present our scenario of critical 
phenomena in gravitational collapse based on the rcnor- 
malization group ideas. To describe our claims clearly, we 
first classify numerically observed (U) universality and 
(S) scaling into two: (1) those of the critical space-time 
and (2) those of the critical exponent. Concerning the 
critical space-time, it has been observed that 

(51) near-critical (but not exactly critical) solutions al- 
ways once approach, but eventually deviate from, 
a self-similar space-time, 

(Ul) which is unique. 

And concerning the critical exponent, it has been ob- 
served that 

(52) For slightly supercritical solutions, the mass of the 
black hole formed is expressed by a scaling law 
Mbh oc (p-pc)''^^, and 

(U2) the value of the critical exponent /3bh is the same 
for all generic one-parameter families in T. 

We could add another kind of universality, universality 
over different systems. It has been found ||l9| that the 
system for some potentials such as F(0) — ^^0^/2, etc., 
also have the above properties with the same /3 as the 
minimally coupled massless scalar field. This suggests 
that the following universality holds: 

(U3a) There are some models which show the same criti- 
cal behaviour (same self-similar solution, same crit- 
ical exponent). In particular, in scalar field col- 
lapse, the critical phenomena is the same for any 
potential V, or at least, for some potentials V. 

Note that the term "universality" is used in two ways. 
In (Ul) and (U2), it refers to the independence on the 
initial data (or a family of initial data), while in (U3a), 
it refers to the independence on the system's EOM (i.e. 
details in the definition of the model considered). In 
statistical mechanics, where the term was used for the 
first time in a similar context, "universality" refers to the 
latter (although there is no "initial data" in equilibrium 
statistical mechanics and thus it is impossible to have a 
counterpart for the former). It is one of our main purpose 
in this paper to give sufficient conditions under which the 
two kinds of universality hold. 



Historically, because of closeness of numerical values 
of the critical exponents of a scalar field a radiation 
fluid 1^, and gravitational waves j^], the following uni- 
versality in still wider sense, though stated less precisely, 
was speculated. 

(U3b) The value of /3bh is the same over many systems, 
including a scalar field, a radiation fluid, and grav- 
itational waves. 

The precise value 0.35580192... in Ref. @ and the val- 
ues for scalar field collapse (0.37) and gravitational wave 
collapse (0.38) shows that this is not true. 

We now present our scenario through two claims, and 
their derivati ons. We do this on t wo different levels. 
First, in Sec. II Bl (resp. Sec. [IB 2| ) we present claims 
concerning qualitative (resp. quantitative) pictures on 
the critical behaviour, and present their (nonrigorous) 
'derivation' based on naive linear perturbation ideas. 
Presentation in this Sec. pB| is intended to provide basic 
ideas, so the following Claims and their 'proofs' are not 
intended to be precise mathematical statements. In par- 
ticular, naive perturbation is not sufficient to deal with 
long time asymptotics rigorously, and this is where the 
"rcnormalization group" idea comes in. We thus present 
in Sec. fl C| more detailed, rigorous derivation of Theo- 
rem |l|, a rigorous version of Claims 1 and 2, based on 
rcnormalization group philosophy. 

A mathematical remark about the phase space is in or- 
der. In order to formulate the time evolution, we consider 
our phase space F to be a real Banach space of functions^ 
C/(^)'s which satisfy certain decay and smoothness prop- 
erties, with a norm || • || under which the phase space is 
complete. Concrete choice of the decay and smoothness 
properties should be determined on physical grounds for 
each models, together with the requirement that F is in- 
variant under the time evolution. The term "approach" 
and "deviate" in the statement (SI) should be interpreted 
in terms of the norm of the phase space, although other 
notion of convergence, e.g. pointwise convergence, could 
equally work. Precise final definition should be given 
based on rigorous mathematical analysis, which is be- 
yond the scope of the present paper. 

In the concrete cases discussed later, the fixed point 
U* does not represent an asymptotically flat space-time. 
Though one is usually interested in asymptotically flat 
space-time in gravitational collapse, we define our phase 
space F to be sufficiently large to include U* . Asymp- 
totically flat data are included in the phase space (if one 



If there are some gauge degrees of freedom, which is often 
the case for general relativistic systems, we should consider 
r as a space of gauge equivalence classes of functions. Or 
we could fix the gauge and consider F to be a usual function 
space. In the latter case, we should be careful in gauge-fixing 
so as not to count "gauge modes" as physically meaningful 
eigenmodes in linear perturbation analysis. 
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chooses the functions Ui appropriately) and form a sub- 
space of r. 

1. Critical space-time 

We first consider the qualitative behaviour of (near-) 
critical solutions. 

Claim 1 Suppose there is a fixed point U* , with no 
marginal modes. If (SI) and (Ul) hold with this U* , 
then there should he a unique relevant mode with a real 
eigenvalue for this U* . Conversely, if U* has a unique 
relevant mode, then (SI) and (Ul) hold with this U* , 
at least for all the initial data close to U* . (Except, of 
course, for the exactly critical initial data.) 

Remark. Because we are considering real PDEs, the 
ODEs 1^ for eigenmodes are also real, and thus all eigen- 
modes with complex eigenvalues (if any) appear in com- 
plex conjugate pairs. This means that if such a system 
has a unique relevant mode, the eigenvalue must be real. 

Remark. (Reason for employing the real formulation.) 
One could proceed without employing the real formula- 
tion. However, in this case, dimension counting becomes 
more complicated; because the desired critical behaviour 
is observed only when the relevant subspace is (real) one- 
dimensional, one has to see whether a given eigenvector 
spans a subspace of (real) dimension 1 or (real) dimen- 
sion 2. By reducing the system of PDEs in its real form, 
we can avoid such subtleties. We emphasize again that 
this does not restrict the applicability of our claims in 
any way. For example, a system of PDEs (for complex 
unknowns) which has exactly one relevant mode with a 
complex eigenvalue is easily seen to have a relevant in- 
variant subspace which is (real) two-dimensional, thus 
leading to the same conclusion (i.e. this does not show 
usual critical behaviour) as is obtained from the real for- 
mulation. Based on this observation, we in this paper 
drop "Re" in the exponent relation /3bh = (3/k. 
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FIG. 1. Schematic view of the global renormalization group 
flow. 

general intersect with the critical surface, and the inter- 
secting point Uc will be driven to U* under the RGT: 

lim \\U^'^ -U*\\= 0. (2.18) 

s — >oo 

Thus, Uc is the initial data with critical parameter pc- 
An initial data Umit in the one-parameter family /, which 
is close to Uc will first be driven towards U* along the 
stable manifold, but eventually be driven away along the 
unstable manifold. Therefore, claims (SI) and (Ul) hold. 

To show the converse, first suppose U* had no relevant 
mode. Then any initial data (sufficiently close to U*) 
would all be driven towards U* until finally be absorbed 
into U*, which contradicts (SI). Next suppose U* had 
more than one relevant modes. Then the stable manifold 
of U* would have codimension more than one, and thus 
a generic one-parameter family of initial data would not 
intersect with the stable manifold of U*. This means 
that a generic critical or near-critical solution would not 
approach the self-similar solution U*; instead, they would 
be driven towards a different self-similar solution with a 
unique relevant mode. Thus claims (SI) and (Ul) do not 
hold. □ 



'Proof. ' Suppose U* has a unique relevant mode. Then 
the renormalization group flow around the fixed point 
U* is contracting except for the direction of the relevant 
mode (Fig. at least in some neighbourhood Af C T 
around U*. Thus there will be a critical surface or a 
stable manifold W^{U*) of the fixed point U*, of codi- 
mension one, whose points will all be driven towards U*. 
And there will be an unstable manifold W"([/*) of di- 
mension one, whose points are all driven away from U*. 
A one-parameter family of initial data I G Af will in 



^PDEs, when one is considering in more than two dimen- 
sions, as remarked at the beginning of Sec. [I A 



Remark. Eq. (2.18) states that the fixed point U* is 
on the closure of the subspace of F which consists of 
asymptotically flat data. This explains, in rather precise 
terms, that the self-similar solution plays such an im- 
portant role that its perturbation gives an exact critical 
exponent even when one is interested only in asymptot- 
ically flat collapses. Asymptotically flat data can con- 
verge to the self-similar fixed point with infinite mass in 
the above sense, even though the ADM mass conserves. 
It should be emphasized that there is no contradiction 
involved here. That is, the mass measured in the unit of 
the scaled coordinate becomes infinite as the space-time 
approaches the self-similar solution, while the physical 
mass measured in the unit of the original coordinate re- 
mains finite. 
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Rema rk. S ome reader might doubt the vahdity of our 



claim ( 2.18 ), because often the critical space-time satis- 
fies Imir^oo o,{r) = aoo > 1- Howe ver, it is not difhcult 
to find norms |1 • ||, for which (2.18) holds; e.g. take 



\a\\ — sup(l + rY 

r>0 



Hr)l 



sup e 

r>0 



\a{r)\. 



(2.19) 



The remaining, highl y non trivial, question is whether the 
Assumptions in Sec. |ll C| are satisfied in this norm, be- 
cause the spectrum of perturbations can depend on the 
choice of the norm (and thus on the function space)Q. In 
other words, the norm and the function space T should 
be so chosen that the Assumptions hold (and that T in- 
cludes U*): if we can find such a norm and F, the critical 
behaviour should be observed; this is what is proved in 
Sec. p C| . The question of the right choice of the norm 
must be answered for each models under proper mathe- 
matical consideration. We emphasize that our goal is to 
present a general framework (physically reasonable suf- 
ficient conditions) under which the critical behaviour is 
bound to happen, which is independent of the details of 
specific models. 

Remark. Fixed points with more than one relevant 
modes are responsible for the so-called multicritical be- 
haviour in statistical mechanics. We give examples of 
self-similar solutions with more than one relevant modes 
(and thus are not responsible for generic critical be- 
haviour) in section VH. Another example is provided by 
the self-similar solutions of complex scalar fields |^^. In 
this model the continuous self-similar solution has three 
relevant modes, and thus is irrelevant for a generic crit- 
ical behaviour. In other words, we have to adjust three 
parameters (instead of one) in the initial condition, to 
observe a critical behaviour governed by this self-similar 
solution. 



2. Critical exponent 

We now turn our attention to the quantitative aspect 
of the critical behaviour, i.e. the critical exponent /3bh- 

Claim 2 // the relevant mode is unique with a real eigen- 
value K, the black hole mass Mbh satisfies 



Mbh (X (p-pc) 



/3bh 



(2.20) 



with 



^Although this dependence, especially on the behaviour as 
r ^ oo, may not be so strong in the concrete case of perfect 
fluid, because the sonic point places an efi'ective boundary 
condition 



stage 




FIG. 2. Schematic view of the flow. Each dot denotes ap- 
proximate location of the flow under iterations of TZa, start- 
ing from several different initial data. Note that the points 
are most densely distributed around the fixed point U* . The 
decomposition into two stages I and II is shown by vertical 
dotted lines . The set A± will be used in our detailed analysis 
in Sec. 



[IC 



/3bh — — 

K 



(2.21) 



for slightly supercritical solutions, where (3 is the scal- 
ing exponent of the scaling transformation (2.i). Here 



f{p) cx (p — Pc)^^^ means there exist p-independent posi- 
tive constants C\,C2 such that Ci < f{p){p — Pc)^^^^ < 
6*2 holds for p sufficiently close to Pc- 



'Proof. ' Th e claim, in particular the exactness of the 
relation (2.21), is one of the main conclusions of this pa- 
per. The derivation of this claim is rather lengthy, al- 
though it is a variation on rather standard treatment 
of similar problems in (rigorous) renormalization group 
analysis of critical behaviour in statistical mechanics. We 
first explain the rough idea, and then explain each step 
in detail. 

We first consider the fate of an initial data U-mit in the 
one-parameter family, which is close io Uc {e — p ~ Pc)'- 



Uc + eF, 



\F\\ 



1. 



(2.22) 



The data will first be driven towards U* along the critical 
surface, but eventually be driven away along the unstable 
manifold, until it finally blows up and forms a black hole, 
or leads to a flat space-time. 

To trace the time evolution, we divide the whole evo- 
lution into two stages: I (s < si) and II (s > si). See 
Fig. ||:^ the stage I corresponds to Umit — > Q, while stage 



^More general situation, where the initial data is not close 
to U* but is close to W(?7*), can be treated in a similar way 
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II to Q ^ (black hole). The stage I is the portion of the 

is'] 

journey where J/jn/t is close to U* and we can use the 
perturbation. Stage II is the region where perturbation 
breaks down, and si is so chosen to make this separation 
possible. Concretely, we choose si as the smallest s such 

that |iC//nit becomes of order 1 (after U^^^^ begins 

to deviate from U*). As will be explained in detail, near- 

critical solution U>^^(^ spends most of its time in stage I, 
and the time of stage II is relatively short. 

Based on the analysis of the flow, we then derive the 
exponent by interpreting the result in our original coordi- 
nate. As will be explained, what matters in determining 
the exponent is the behaviour of the flow in stage I. 

Step 1. Flows in stage I: 

The analysis of this stage is essential, because the flow 
spends most of its time at this stage, and gives the value 
of the exponent. 

For Stage I, we can rely on linear perturbations to get 



e e 



KSi 



F, 



rcl 



(2.23) 



where i^roi is roughly the component of the relevant mode 
in F. The derivation of (2.23) requires some calculations, 
together with suitable assumptions on invariant mani- 
folds of T. For the ease of reading, we postpon e thi s 
derivation to section [IC3, and proceed assuming ( 2.23| ) 
is correct. 



Due to (2.18), we have (for large si) 



U 



^ F 



Our choice of ,si implies 



(2.24) 



(2.25) 



(Here 1 means a positive small which is independent of 
e.) Note that s\ diverges logarithmically in e as e ^ 0. 

Step 2. Flows in stage II: 

We further evolve U-^^^'^ until it finally diverges or leads 
to a flat space-time. At this point, we assume that the 
phase space F has a rather simple structure, i.e. all points 
(at least in a neighbourhood oft/*) on one side ofW^{U*) 
blow up and lead to a black hole, while those on the other 
side disperse and lead to a flat space-time. Once one 
admits this nontrivial assumption on the global structure 
of F, stage II becomes relatively trivial, because of the 



large second term of (|2.24|). That is, the data U^^l^ differs 



from U* so much that one can tell the fate of this data 



as long as the flow along ^{U*) is regular. This is because 
the time spent before it reaches Stage I is finite. 



depending on the sign of e. In particular, if a black hole is 
formed, the solution will blow up in a finite, rather short, 
time (uniform in e) . This means the radius of its apparent 
horizon, and thus its mass, will be 0(1) measured in ^. 

Step 3. Going back to the original scale: 

Finally, we translate the above result back into our 
original coordinate (i, r). The relation r = ^e"^" implies 
that the radius of the apparent horizon, which is 0(1) 
measured i n i s in fact 0(6"^"') measured in r. So we 
have from ( 2.25| ) 



Ah 



BH 



0(e-'3^') = 0{\ef'^) 



(2.26) 



Therefore the critical exponent is given exactly by 



/3 



(2.27) 



This concludes our 'derivation.' □ 
Remark. As can be seen from the above 'proof,' an 
intuitive understanding of the exactness of /3bh is as fol- 
lows. In the limit e — s- 0, the time spent in a finite region 
around U* diverges whereas the time spent outside is 
usually finite before a black hole is formed, if there are 
no other fixed points, etc near the fiow. So the change 
of the profiles of the fields are dominated by the linear 
stage so that the exact critical exponent is given in terms 
of the eigenvalue of linear perturbation aro und U* . The 
rigorous argument is presented in Sec. [IC. 

It follows immediately from the above claim that con- 
ditions (S2) and (U 2) hol d and the critical exponent /3bh 
is exactly given by (2.27). The above 'proofs' show that 
our claims are true not only in asymptotically flat cases 
but also in cases with any definition of mass proportional 
to the typical scale of the functions C/(C). 

If one admits Claim 2, it would in fact be a big surprise 
if the condition (U3b) did hold, since the critical solution 
of other systems, such as that of a scalar field system 
and that of a gravitational wave, are so different that 
one cannot expect the eigenvalues of the perturbations 
around them to be the same. An explicit example of 
nonuniversality is given by the critical exponent /3bh for 
various par ameter values of 7 for the perfect fluid: see 
section VII. 



C. Detailed analysis of renormalization group flows 

Simple linear perturbations are not sufficient for rigor- 
ous study of the long time asymptotics of the solution. 
In this section, we describe the "renormalization group" 
type analysis in detail, which goes beyond the linear per- 
turbation and provides firm grounds on our claims. Our 
goal is to derive Theorem which is a rigorous version 
of Claims 1 and 2, from a few "assumptions." The "as- 
sumptions" from which we start are grouped into two: 
(1) assumptions on the linear time evolution of the sys- 
tem, defined by T (or equivalently 7^,t/» — e"'^), and (2) 
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certain regularity (smoothness) properties of TZs and Tg 
which physicaUy look innocent. 

In the following cr is a positive constant suitably chosen 
for both of the following assumptions to hold. The first 
assumption specifies the properties of the linear part of 
the time evolution, generated by To,u*'- 

Assumption LI (Invariant subspaces of Xj^u*) For 
a suitably chosen a, the whole tangent space Tu*T is a 
direct sum of invariant subspaces of T^^jy* '■ 



(2.28) 



Here E^{U*) and E^{U*) are invariant subspaces of 
7^ [/• , of dimension N and of codimension N, respec- 
tively. There exists k > such that the restrictions of 
%j^U' on these invariant subspaces satisfy: 



\\T,M4F)\\>e'^-\\F\\ 
\\%,U'iF)\\<e---'^\\F\\ 



(Fe £;"([/*)) 

(F ££;"([/*)) (2.29) 



holds with some (F-independent) Kq > 0. The first in- 
equality of d^H) is due to completeness, while the second 
is nothing but the triangle inequality. Because this box 
norm simplifies our presentation (e.g. conditions on var- 
ious constants in Propositions), we extensively use this 
box norm in the following, writing it simply as || • || . Some- 
times, we write \\F\\s = ||Fi, = 

Notation We exten sively use following abbreviations 
throughout Sec. |lC| to Sec. pE|. 



• We write a V 6 = max{a, b}. 

• 0(a) denotes a scalar or vector depending on the 
context, whose magnitude (measured in absolute 
value I • I for a scalar, in the box norm || • || for a 
vector) is bounded by \a\. 

• We also write f ^gzthto simultaneously represent 
f < g + h and / > g_— h. That is, / > g ± /i is 
equivalent to f = g + 0{h) when h > 0. 



The estimates ( 2.29| ) in Assumption LI could be fur- 
ther deduced from suitable assumptions on the spectrum 
of 7(7*, because Ta-^jj* = exp(s7(7. ). One possibility of 
such an assumption would be: 

The spectrum^ of T = 7(7* is a union I]_ U 
E+. Here (i) S+ is a finite set with total mul- 
tiphcity TV, (ii) S_ C {z £ C | Rez < -R'} 
and S+ C {z G C I Rez > k'} for some k' > 0, 
and (iii) T is a sectorial operator. 

For a related estimate, see e.g. Henry Theo- 
rem 1.5.3]. 

Before stating the second assumption, we set up the 
following convention on the choice of the norm || • ||. 

Convention on the choice of the norm | • | . For a 

vector F, we denote by F^ (resp. F") the E^{U*) (resp. 
F"(C/*)) component of F: 



F = F"" + F'' 



(2.30) 



with F"" e F"([/*), F" e F"(J7*). Assumption LI guar- 
antees the possibility of the decomposition. We define 
the "box norm" 



|F||box^sup{|lF«|l,|lF"|l}, 



(2.31) 



Note that this box norm is equivalent to the original || • ||, 
because 

Ko\\F<^lKo\\F-\\<\\F\\<\\F^\\ + \\F-\\ (2.32) 



The second assumption refers to regularity properties 
of the flow, and essentially requires that the map TZcr 
has a Lipschitz continuous Frechet derivative. Although 
it is sufficient to require only Holder continuity (of the 
derivative) with a positive exponent in the sequel, We 
here require Lipschitz continuity, to make the presenta- 
tion as simple as possible. 

Assumption L2 (Smoothness of TZa) There exist pos- 
itive constants 6q, Ki, K2 such that 

(i) the formal relation (|^) holds for || [/-[/* ||, ||F|| < Jq: 

naiU + F)=n„{U)+%^uiF) + 0{Ki\\F\\^) (2.33) 

(ii) Xj^jj is uniformly bounded for ||?7 — < Sq: 

I17;,C7(F)|| <if2||i^||. (2.34) 
We take 60 sufficiently small so that SqKi < K2- 

Note that Assumption L2 implies the following 
Lipschitz- type continuity for TZ^ and %ju (for \\Ui — 
U*\\<So {i = l,2)): 

n^iUi) - 7^,(C/2) = 0(2X2 \\Ui - C/2II), (2.35) 
(r,,(7i - %^u,) F = 0{QKi -WUi- C/2II • ||F||). (2.36) 



Here we of course mean both continuous and discrete 
spectrum. 



[To prove ( 2.3(]| ) from Assumption L2, first note it suffices 
to consider the case of ||F|| = \\Ui — C/2II/2, because T^^u 
is a linear operator. Then evaluate Ai? = TZ.^{U2 + F) — 

nAU2)'nAUi + F) + nAUi) for ||F|j = - [/2II/2 

in two ways using (2.33): (1) express it in terms of 7^^(/;, 
(2) expand each Tlcr{-) at {Ui + U2)/2 and show that ATZ 
isO(10Xi||F||2).] 

Although Assumption L2 physically looks rather inno- 
cent, its verification is a nontrivial mathematical problem 
depending on individual cases. We emphasize that Ki 
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term in Assumption L2 is not necessarily bounded uni- 
formly in cr, and can diverge for large a. For example, 
they will certainly diverge (as cr — s- oo) for non-critical ini- 
tial data which will lead to a black hole. What is claimed 
in the Assumption is that one can find finite nonzero a 
and (5o for which the above claims hold. 

We follow the time evolution by considering a discrete 
dynamical system generated by TZna {n G Z,n > 0). 
The reason for considering the discrete system is twofold. 
First, it is naturally generalized to the case with discrete 
self-similarity, as discussed in Sec. [ID. Second, the op- 
erator TZ^ could be more regular than TZ, thanks to a 
smoothing effect of the time evolution over a nonzero pe- 
riod cr. (I.e. the linearization T may not be bounded, 
while T can be.) 



1. Local behaviour of the flow: Existence of the critical 
surface (stable manifold) 

Based on the above assumptions, we now start our 
renormalization group study of flows. We here consider 
the qualitative aspects of the flow in a vicinity of the 
fixed point U* . 

We remark that there are considerable number of 
mathematical literature concerningthe local structure of 
invariant manifolds. Proposition |l| is a celebrated "in- 
variant manifold theorem"; see e.g. Shub Theo- 
rem 5.2, Theorem II. 4] or [g2[ Theorem III. 8], for precise 
statements and a proof by sophisticated methods under 
slightly different assumptions. Although the cited theo- 
rems are sufficient for our needs in t his s ection, they do 
not cover the case considered in Sec. 



[I E , where we con- 



sider the universality over different models. We therefore 
present our version of an elementary proof, in order to 
make the presentation more self-contained, an d als o to 
pave a natural way to our later analysis in Sec. p^I E| . 

We consider the time evolution of an initial data U 
(close to U*) under TZna, by expressing it as {n > 0) 



f/("")=7^„,(C/) = C/* + F„ 



(2.37) 



where the second equality defines the difference Fn be- 
tween TZncriU) and U*. See Fig. ||(a). Assumption LI 
guarantees that we can always decompose this way. 
We begin with the following lemma, which claims that 



an initial data with lli^nl! < 



is driven away from U* 



(at least until becomes somewhat large). 

Lemma 1 Define Si > by 



Si = min < (5o 



1 - e- 



2Ki 



(2.38) 



and consider the time evolution of an initial data U = 
U* + Fq with ||i^o||s < ll^ollu < 5i. Then as long as 
||-Pra||u < 5i, we have 










/ ^ (/"" 




1 ^ 

/ " ^ 




£" 




(b) 



FIG. 3. Schem atic v iew of three decompositions: (a) the 
decomposition of (2.37), where F measures the difference be- 
tween U and [/*, and is deco mpose d resp ect to and (b) 
the decomposition of (2.77), and (2.78), where Fn measures 



the difference between [/f"""' and and then is decom- 

posed with respect to E'^ and (c) the decomposition of 
(|A7|), which decomposes with respect to and VK". 



\\Fnh < \\FnU (2.39) 
e''"/' < ||F„+i||u < 2X2||F„||u (2.40) 

Proof. We proceed by induction in n. We start by 
obtaining recursion relations for Ff^ and F". Assump- 
tion L2 reads 

n^{U* + F„) = n,{U*) + %M'iFn) + 0{Ki\\Fj^) 

(2.41) 

which implies (note: TZa{U*) = U*) 

Fn+i ^ %,U'{Fn) + OiKiWFj^) 

-%^U'{F;^) + OiKi\\Fr) 



Xt,U' {Fn 
'Fa,U' {F^ 



(2.42) 



where in the second step we used the linearity of Ta-^u* 
and the fact that F^ and F^ bel ong to invariant subspaces 
of Xj,u'- Decomposing as in ( 2.37 ), we have recursion 
relations for F,' (i = s or u) : 



Fn+l 



%,u*iF^) + 0{Ki\\Fn 



(2.43) 



Our inductive assumption ||F^|| < ||f"| | imp lies ||F|| = 
||F"|| , by the defin ition of the box norm (2.31). So taking 
the norm of dHI) and utilizing Assumption LI and As- 
sumption L2 to bound ||7^ .[/* (F^)||, we have recursions: 



F„V||<X2||F„"||+Xi||^^,^||2 



(2.44) 
(2.45) 
(2.46) 



Now the recursion (2.44) immediately implies (for 
!i^„"|| < <5i < ^0 < K2/KT) 

iif;viIi < {K2+Ki\\F::\\)\\F:i\\ < 2k,\\f;:\\. (2.47) 



On the other hand, ( |2.45| ) implies (for < Si] use 

also the definition of 5i) 
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(2.48) 



These two prove ( |2.40| ). [Note: (1 + a)/2 > ^ for posi- 
tive a.] 

Finally, ( ^.46 ) implies 



n+ll 



<e-^'^\\F^:\\+K^\\F,':r 



(2.49) 



This, combined with < just proven, implies 



( ^.39 ), because 



KiSi < 1. 



Proposition 1 Under Assumptions 1 and 2, we have 
(i) [Existence of critical solutions.] For any F'^ G E^{U*) 
with WF^'W < 6i, there exists F^"(F") e £:"([/*) such that 



Moreover, 



F^iF"") -^U* as n 



(2.50) 



\f:{f^)\\<c^\\f 



IS II 2 



with Ci = 



2Ki 



1 - e- 



(2.51) 



We define the "local stable manifold" as 

W^,JU*) = {U* + F' + F^iF^)\\\F^ <S,}. (2.52) 
(ii) [Motion on Wf^^{U*).J There exists k> such that 
WR-naiU) -U*\\< e-"^'"||C/ - U*\\ (2.53) 

and 

\\n(^n+iy{U) -U*\\< e-^'^ |l7^„,(C/) - U*\\ (2.54) 

for any vector U € W^ioc(^*) which satisfies \\U — U*\\ < 
di. We can in fact take R ^ R by taking \\U — U*\\ small. 
(Hi) [Invariant manifolds.] Above W[^^(U*) in fact forms 
a manifold of codimension N . 

We only present a proof for (i) , (ii) of the Proposition. 
Detailed information presented in (iii) is not necessary 
for our later analysis, so we only make brief comments 
on their proof at the end of the proof of (i) and (ii). 
Our proof is modeled after the argument by Bleher and 
Sinai | p3| , which first appeared in a similar context in 
statistical mechanics. 

Proof of Proposition [J, (i) and (ii) 

We mainly consider the case of iV = 1, which is of our 
main interest. For the other cases, see the end of this 
proof. We denote the normalized relevant mode of T„_u' 
by i^"': %,u*F'"'^ = e'^'^i^"', = 1. 

The proof is done in several steps. 

Step 1. The goal. We are interested in the time 
evolution of vectors U^^^ = U* + F^ + aoF"'\ for fixed 
Fq and various ag. Our goal is to show that there is at 



least one value of oq (called function of Fq) for 

which lim„^oo7^„„(C/(°)) = U*. 

Step 2. Recursion relations for a„ and F^. We ob- 
tained recursi ons f or F^ and F^ = a„F"' in the proof 
of Lemma 1^: ( |2.43 ). Utilizing Assumption LI, and using 



|F|| < max||a„|, I F^\ \} which follows fr om th e definition 
of the box norm ( |2.3l| ), we can convert (2.43) into recur- 
sion relations for a„ and /i„ = \\F^\\. The result is (recall: 
ay b = max{a, 6}) 



a„+i = e'"^a„ + 0(^1 (a„ V /i„) ) 
K+i < e-'^^hn + 6{Ki{an V /i„)2) 



(2.55) 
(2.56) 



Roughly speaking, the recursion shows that the map 
for a is expanding, while that for h is contracting, with 
"shifts" of 0{Ki{an V /i„)^). If we neglect the sec- 
ond term, a„ would blow up. Thus the stable set 
W[^^{U*) is located by carefully balancing the inhomo- 
geneous 0{Ki{an V /i„)^) term with the first expanding 
term, so that a„ will remain bounded (in fact goes to 
zero) as we iterate the recursion. That we can in fact 
find at least one such ao(= flc) for given ho can be shown 
by using the argument due to Bleher and Sinai as 
explained in Steps 3 — 5. 

Ste p 3. Continuity of a„ and hn as functions of Oq- 
from (2.35) TZ^ is a continuous function of its arguments. 



This implies in particular that a„+i and F^^^ are contin- 
uous functions of a„ and F^. Iterating this for n times, 
we see that a„ and F^ are continuous functions of the 
initial data, qq, for fixed Fq. This implies hn = \\F^\\ is 
also continuous in oq. 

Step 4- Solving the recursions. We now solve the 
recursion, following the technique of Bleher and Sinai. 
We will show that there are positive sequences r„, t„ such 
that 



as n 



oo 



(b) < hn^i < rn+1 holds, as long as |a„| < i„ and 

< hn < rn. 

(c) As long as < /i„ < r„, the set of all a„+i when a„ 

sweeps the interval (— tn,tn) contains the interval 
(-t„+i,t„+i) inside. 

The proof of (a) — (c) is postponed to "Step 5" for the 
ease of reading. 

Now given (a) — (c), we can show the existence of the 
critical Oc as follows. Let /„ be the interval /„ = {z £ E : 
I z I < tn}. Because a„ is a continuous function of oq , we 
can consider the inverse image Inv„(/„) of the interval 
/„ under the n-iterations of the recursion. But (a) — (c) 
show that Inv„(/„) are nonincreasing (in n) sets, with 
nonempty intersections. So, taking an initial condition 
Q-c € n„Inv(/„), we can guarantee for all n > 



< tn and hn < rn. 



(2.57) 



Because r„ , i„ ^ as n oo , ( |I57| ) shows U* 
for this special a^, i.e. U* + F"" + F^ G W^iU*), with 
F" ~ flcF™': existence of a critical solution for given F^. 
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Step 5. Proof of (a) — (c). Now we give an example 
of sequences r„,<„ which satisfy (a) — (c) of Step 4. We 
define r„, i„ according to (for < rg < 5i) 



tr. 



Cirl 



with Ci 



2Ki 



(2.58) 
(2.59) 



Then (a) is easily seen to be satisfied, because e '"^ + 
KiTo < (1 + e-^'^)/2 < 1 by the choice of 5i. To prove 
(b) and (c), first note the following sufficient condition: 



r„+i > e ^'"r„ + ifi(r„ Vi„) 



(2.60) 
(2.61) 



The choice of Ci guarantees Ci^n < 1, and thus t„ < 
rn- T his, t ogeth er with k > R, immediately shows that 
( b.60 p and ( ^.6lD are satisfied. Eq.(|2.59D for n = proves 



S'iep Motion of the critical solution. Now that 
we have proven the existence of a critical solution, we 
concentrate on its time evolution. We first note |a„| < 
hn holds for the critical solution for all n > 0, because 
otherwise the solution is shown to be dr iven away from 
U* by Lemma |l|. But then, the recursion ( ^.56| ) is reduced 
to 

K+i<[e-'^'' + KiK]K (2.62) 
which, in view of hn < r^i < rQ, implies 



\n^n+i)Au)-u*\\ < e-^^ ||7^„„ ([/)-[/* II 



with 



Kiro 



< 



< 1 



(2.63) 



□ 



Cases of N ^ I and a comment on Proposition ( Hi ) 
All the above proves the existence of W^^^{U*) and the 
estimate (2.54) for the case of = 1. For N > 2 we can 



proceed similarly and find the critical surface (this time, 
we have to adjust N parameters). For A'^ = 0, similar 
argument shows that the map TZ^ is a contraction in a 
neighbourhood of U*, thus proving that all the initial 
data close to U* is attracted to U* . 

Proof of Proposition |l|, (iii) can be performed, e.g., as 
in PI, Theorem II. 4]. □ 



2. Local behaviour of the flow: Unstable manifold 

In the previous section, we established the existence 
of the critical surface (stable manifold). We also showed 
that initial data which is sufficiently away from the crit- 
ical surface (characterized by \\U - U*\\s < \\U - U*\\u) 
is driven away from U*. In this and the following sec- 
tion, we concentrate on the fate of near (but not exactly) 



critical solutions. We in thi s sect ion define an "unsta- 
ble manifold." Then in Sec. II C 3, we show that all the 



near-critical initial data (starting from an neighbourhood 
of [/*) are at first driven towards U* along the critical 
surface, but are eventually driven away along an unstable 
direction. We also give quantitative estimates on time in- 
tervals until the deviation of the data from U* becomes 
large. 

We concentrate on the case where there is a unique 
relevant mode (A^ = 1 in Assumption LI). More precisely 
we assum e, in addition to Assumptions LI and L2 of 
SecllI C 1|, 



Assumption LIA (Uniqueness of the relevant 
mode) The relevant mode of Xj^w is unique, i.e. A = 1 
in Assumption LI. We denote the relevant mode by 
and the relevant real eigenvalue by k (> k): 



(2.64) 



Throughout this section, we use following notation 
• 62, (^2 are positive constants chosen so that 

1 - 



62, £2 < min < Si 



S2 + £2 < 



\2Ki 



l2Ki 



(2.65) 



(2.66) 



• Given U GT, define ni(J7;e) to be the largest in- 
teger such that \\TZkcr{U) ~ U*\\ < e holds for all 
< A: < ni{U;e). I.e. 

7ii(C/;e) = max{r^|||7^fc^(C/)-C/*|| < e for k < n} 

(2.67) 

To trace the time evolution, we have to specify the 
unstable manifold. If TZ^ is invertible, we can define the 
unstable manifold as 



W{U*) = {U\TZna{U) ^ U* 



)} , (2.68) 



interpreting TZna — (TZ-na)^^ for negative n. However, 
we here do not assume the invertibility of TZ^. So the 
definition of W^{U*) requires some additional work. Al- 
though our definition is somewhat artificial, it suffices 
our need. 

We define several sets, including the local unstable 
manifold W^^{U*). Although all of these are defined 
with respect to the fixed point U*, we do not explicitly 
write U* except for VF["^(C/*), to simplify notation. We 
first introduce a line segment Dr for r e R as 

^ {(1 - t){U*+rF'''^) + tTZa{U*+rF'''^)\0 < t < 1} 

(2.69) 
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and define Wr as [ni is defined in ( ^!67| )] 

Wr = {nnaiU)\U eDr,n< m ([/* + rF'^l; £3) } . 

(2.70) 

Then we define the local unstable manifold at U* as 



W,UU*) = lim 

r\0 



WrUW- 



(2.71) 



Lemma ^ below guarantees that W^^(U*) is well de- 
fined, and all flows close to U* are squeezed around this 



Lemma 2 (i) [Structure of W^^^iU*)] W^^^iU*) of 
(2.71) is a well-defined set, and is "one- dimensional" in 
the sense that its intersection with a plane U* -\-E^{U*)-\- 
dprci exists as a unique point for \a\ < £2- Moreover, any 
U,U'eWiUU*) satisfy 



\U-U'\\s < Ci\\U-U'\U 



with 



Ci = 



7X162 



- 7i^ie2 



<1 



(2.72) 



(2.73) 



(a) [Contraction around W^^{U*) inE^-planeJ Consider 
the time evolution of a vector U under TZna, with the 
initial data satisfying \\U—U*\\ < £2- We can find R" > 
(depending on £2) such that 

dist(nnAU),W,UU*)) < e-"«"n|[/- [7*11 (2.74) 

holds as long as \\TZna-iU) — U*\\ < £2- 

For a proof, see Appendix ^. 
Finally we define sets of vectors 



A. 



e r 



-2k<7 



< \\U±~U*\\,,<e, 



dist([/±,T4^i^',(C/*)) <<5}, (2.75) 

where for a given set A C T and a vector U E T, we 
denote their "distance" as 



dist(J7,yl) = inf \\U-U'\\ 
U'eA" 



(2.76) 



The set A± is schematically shown in Fig. ^of Sec. II B 2 



3. Local behaviour of the flow: Near-critical solution 



the essential role in determining the critical exponents. 
Our presentation is modeled after a standard approach in 
applications of rigorous renormalization group techniques 
in statistical mechanics, see e.g. [P4lP5|. The rest of the 



ec. |II C 4 



flow (Stage II) is considered in Sec. |II C 4 . 

Time evolution of an initial data sensitively depends 
on its location relative to Wl^^{U*). So we decompose 
the initial data U as 



U = bF' 



rcl 



u. 



(2.77) 



where Uc € W{^^{U*) [i.e. 7^„<,(^7c) ^ J7* as n ^ 00], 
F'^'^^ is the normalized relevant eigenmode deflned in 
( ^.64, an d & S K is its coefficient. See Fig. |(b) in 
Sec. lie 1. Proposition |l| (i) guarantees that we can find 
Uc and decompose this way for \\U — U*\\ < Si. Generi- 
cally b — 0(p — pc), because b measures the component 
of mU-Uc. 

Now, for n > 0, we define F„ in terms of the decom- 
position ^ 



nna{U)=nnaiU,)+F„ 



and further decompose Fn in the spirit of ( 2.3Cl| ) 



F„, 



br,.F 



rcl 



F^ 



(2.78) 



(2.79) 



See Fig. ||(b) in S ec. II CI . We do emphasize that 
the decom position (2.78) differs from that of (2.37) in 
Sec. II C 1 , in that F there measures the difference be- 
tween U and the fixed point U*, while F here measures 
the difference between U and the (time-evolved) critical 
solution C/^'"^^. Also note that F„ (n > 0) is not neces- 
sarily parallel to F'^°^. 

The initial condition implies bo — b, Fq — 0. 

After these preparations, concerning the flow in Stage 
I, we have: 

Proposition 2 (Local behaviour of non-critical flows) 

Consider the time evolution of an initial data Uinit = 
6oF>-<=i -f C/c, with U, £ W^^^iU*), ||C/c - U*\\ < 62, 



< Ibol < £2; o-nd at each na we decompose as in ( 2. 7^ : 
7^„o■(C/i„it) = TlnaiUc) + 6„F'<=i -f F,l. Let = ma^TTe 
2| \bk\ < £2 for < fc < n}. Then under Assump- 
tions LI, L2, and LlA, we can find C2 > 1 (independent 
of bo and €2 ) such that the following holds: 
(i) n2 is finite and satisfies (n < n2) 



3e«'^ - 1 
F^A\ < \b„ 



e2< |6„J <e2, 



(2.80) 
(2.81) 



We now trace the time evolution of a slightly super- 
critical initial data U. Our final goal is to trace it until it 
finally blows up and forms a black hole, but in this section 
we trace it until its deviation from Wi^JU*) becomes of 
order one (Stage I of Sec. II B 2| ) as in (2.24), so that we 
can tell that it is not critical. As shown, this part plays 



^ We could trace the time evolution by employing the de- 
composition ( 2.77 ) also for n > 0, i.e. TZnaiJJ) = Uc n + b nF''°'^ 
with Uc,n £ W^iocC^*)- However, we here employ ( 2.78), b e- 
cause it also allows us to deal with the case of Sec. [IE in 
essentially the same way. 
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||7^„.(c/< 

£2 



C-2 



U*\\ <e2 



< C2e2. 



1 + e« 



(2.82) 
(2.83) 



(ii) For sufficiently small \bo\, we can find S > such 
that 



(2.84) 



where plus or minus sign is chosen according to the sign 
of bo- In fact we can take (5 \ as |5o| \ 0. 

The above Proposition gives an estimate on the time 
interval spent by the flow in Stage I. It also provides 
information on the exit of Stage I, by showing that the 
exit is essentially squeezed into the tiny tube^ A± around 
W{UU*) (see Fig. |). 

Proof of Proposition (i) We first derive recursion equa- 
tions for bn and 

Using ( p. 33 ) with U = Uc, we have 



7^(„+l),(^7) = 7^,([/("'^)) = 7^,(c/("") + i^„) 
= 7^,(^7("'^)) + r (F„) + o(i^i||^^„in. 



(2.85) 



This implies, according to our decomposition (2.78), 
F„+i = r ^^<„„, (F,,) + 0{K4F^\\^) 

= W(F«) +0(6i^i||t/i"'^) - U*\\ ||F„||) 
+ 6{Ki\\F„f) 



(2.86) 



where in the second step we used ( p^ ). This corre- 
sponds to the first line of ( 2.42| ); the difference is we here 
have 0(6 JCiH [/c - U*\\ \\F\\). Arguing as we did in deriv- 
ing ( |2.43|) — ( ^.46 ), i.e. decomposing F into F^ and i^" 
as in ( p. 79 ), we this time obtain 



fn+l < e fn 

+Kii\b^\ V /„) [0(65„) -I- 0(|6„| V /„)] (2.87) 
bn+i = e'^'bn 



-KiQbr^l V /„) [0(6.g„) + 0(|6„| V /„)] . (2.^ 



|C/< 



(ncr) 



-U*\\,fn=\\F^ 



where we have introduced g. 
to simplify the notation. 

Our remaining task is to solve the recursions. This is 
an elementary exercise in calculus, and is summari zed a s 
Lemr na pl in Appendix [bI the r ecur sive in equ alities (2.87) 
and ( ^.8"^ ) are of the form of (BS) and (B4) there, with 



A 



A 



and C3 = 6K1 . And that gn satisfies 



t he as sumption of the Lemma follows from Proposition]^, 
(^.53|) . Part (i) of Proposition ^ thus follows immediately 
from Lemma a. 

(ii) We apply Lemma || to the time evolution of the 
vector U2 = f7("20-/2)^ Here n2/2 should be interpreted as 
its integer part. We have, by Lemma ^ and Proposition [|, 



< £2 



-(n2-n) 



526 



(2.89) 



For n G [712/2,712], the above is bounded uniformly 
by 62 + (2- This also shows that the initial data U2 
for this part of time evolution satisfies ||t/2 ~ = 
|j[/(»2<T/2) _ [/*|| \ as ri2 / 00. So taking 62 + £2 suf- 
ficiently small, we can apply Lemma to conclude that 
dist(C/'^"2'^\ W{^^{U*))goes to zero as 722 goes to infinity 
(i.e. bo goes to zero) Q □ 



4. Global behaviour and physical consequences 

We are now at the stage of presenting the rigorous ver- 
sion of Claims 1 and 2, together with necessary assump- 
tions. After the analysis in the previous section, that 
is, the analysis of the flow at Stage I, we in this section 
consider Stage II of the flow. We want to trace the flow 
further until it blows up, and then calculate the mass 
of the Black hole. This stage involves nontrivial physi- 
cal and mathematical processes whose detailed analysis is 
beyond the scope of this paper. I n this sense, this section 
is essentially a repetition of Sec. II B 2 , with a difference 
that we here explicitly list up all the assumptions (which 
are physically reasonable in view of numerical simula- 
tions) on the global behaviour of the flow, under which 
we can present the rigorous version of Claims 1 and 2, 
i.e. Theorem |^. 



Assumption Gl (Flows in Stage II) For s ufficiently 
small (^3 > and for A±{e2, S3) defined in ( 2.75| ) we have: 



u'^^ = TZs{U-) tends to a flat space-time as s 
00, for all U- e A_(e2,<^3), 



is) 



TZs{Uj^) blows up at some flnite time s 



SAH > 0, ^ = ^ah; sah and $ah arc uniformly 
bounded above and below for all U+ E ^+(€2, (^3). 



[ Rec all that 1^ is related to 

(PD 



by r = 



-13s 



as m 



*With an additional assumption on the smoothness on Ts^u 
for s G (0, a), we can reduce the tube A± to a disk, given by 

{u± er\ ||[/±-[/*||, = e2,dist((7±,w"(f/*))<5}. 



^Precisely speaking, |6„2l < £2 does not guarantee ||(7'"^'^^ — 
C^*||u < £2, because 6„2 measures [[C/'"^'^' — I7i"^'''||, and in 
general ||f/i"^'^'' — t/'*!! > 0. However, because \\Uc"^'^^ — 
U*\\ goes to zero as 722 goes to infinity, the difference can be 
neglected as bo goes to zero. 
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Assumption G2 (Mass inside the apparent hori- 
zon) The mass -Mah inside the apparent horizon at its 
formation is related with its radius tah by C^^tah < 
Mah < C^AHj where C is a finite positive constant in- 
dependent of TAH- 

Assumption G3 (Relation between masses) The 

final mass of the black hole A/bh, and the mass at the 
formation of the apparent horizon Mah, are of the same 
order. More precisely, there exists a positive constant C" 
such that C'-^Mah < Mbh < C'Mar holds for all the 
near-critical solutions. 

Remarks on the meaning of these assumptions are in 
order. 

Remark. 

1. Assumption Gl specifies the behaviour of the flow 
in Stage II. Thanks to our analysis in Stage I, 
we have only to know the fate of the fiow which 
is emerging from a very sma ll neig hbourhood of 
WiU*) (recaU Fig. | of Sec. |IIB2D . This means 
that physically it is sufficient to trace the time evo- 
lution of two vectors which is located on W^{U*). 

2. Assumption G2 is valid for a spherically symmet- 
ric gravitational collapse (in fact, Mah = ?'ah/2), 
which is the main topic in this paper. However we 
termed it as an Assumption, in order to make our 
presentation more general. 

3. Assumption G3 involves a subtle physical process 
after the formation of an apparent horizon. Strictly 
speaking, the radius of the apparent horizon is not 
exactly related with the final Bondi-Sachs mass 
which gives the mass of the black hole. The fi- 
nal mass Mbh differs from (greater than, when the 
dominant energy condition is satisfied,) the mass at 
the apparent horizon Mah by the matter absorbed 
into the black hole after the formation of the ap- 
parent horizon. However, in our case, it is expected 
from the form of the relevant mode and the blow-up 
profile in numerical simulations that this effect can 
only make the black hole slightly heavier, say by a 
factor of two or so. (It should be noted that usual 
numerical simulation assumes Assumption G3, and 
identifies the radius of apparent horizon with the 
final mass of the black hole.) 

Under these preparations, we can now state our rigor- 
ous version of Claim 1 and Claim 2, which is the main 
result of this section. 

Theorem 1 Under Assumptions LI, L2, LlA and As- 
sumptions Gl, G2, 03, the universality (Ul, U2) and 
scaling (SI, S2) of Sec. hold. More precisely, 
(i) All the initial data sufficiently close to U* (but not 
exactly on Wf^^{U*)) once approach, but finally deviate 



from, U* . The final fate of these data are either a flat 
space-time or a black hole, depending on which side of 
^ioc(^*) initial data. 

(a) Consider the time evolution of the initial data 



bF 



rcl 



(2.90) 



with \\Uc — t^*|| < S2, < |6| < 62, which is on the side of 
black hole ofWf^^{U*). Then as b goes to zero, the mass 
of the black hole formed satisfies 



C,\bf/^<^<C2\bf/- 

Mq 



(2.91) 



with b-independent positive constants Ci, C2. Here Mq is 
the initial (say, ADM or Bondi-Sachs) mass. |^ 



Remark. By refining the above argument leading to 
Theorem |l|, especially refining the proof of Lemma ^ we 
can see the following: if there is a marginal mode, but 
we still have a stable manifold of U* of codimension one, 
it can happen that the scaling of mass has a logarithmic 
correction term, like 



Mbh (X \p* - pI''^" (In \p* - p\ 



(2.92) 



with some exponents /3bh and /Jjn. 

Proof. Proposition H is a ri gorous version of Step 1 of the 
'proof of Claim 2 of Sec. pB^ . Assumptions Gl, G2 
and G3 enable us to perform Step 2, proving (i) of the 
Theorem. This also shows 



rAu/C < Mbh < CrAH, £2/^ < |5|e"" < Ce 



2, 
(2.93) 

where C is a positive constant, and t is the total time 
until the formation of the apparent horizon. Since 
tah/Mo = 0{e^'^) by the definition of scaling transfor- 
mation (p^), we get (ii) of the Theorem. □ 



D. Discrete self-similarity 

Significance of discretely self-similar solutions can be 
well understood with minor modifications to the scenario 



presented in Sec. [I B 



^'^We have included Mo in the expression (2.91) to avoid a 
nongeneric situation in which Mo behaves like power of 60 and 
/3bh becomes different from /3/k. 
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FIG. 4. Schematic view of fiows and the Poincare section 



dist([/,f/ie) =inf lit/ II, (2.96) 



and we simply say "[/ is close to Uic" when dist(J7, Uic) 
is small. 

We can now show the discrete versions of the claims in 



Sec. II B 



Claim 1' Suppose there is a limit cycle Uic of periodic- 
ity A, with no nontrivial marginal modes. If (SI) and 
(Ul) hold with this limit cycle U\c, then there should be a 
unique relevant mode for this U\c. Conversely, if Uic has 
a unique relevant mode, then (SI) and (Ul) hold with 
this Uic, cit least for all the initial data sufficiently close 
to Uic- (Except, of course, for the exactly critical initial 
data.) 

Claim 2' // the relevant mode is unique, the black hole 
mass satisfies 



A/bhcx(p-p*)'3/- 



(2.97) 



1. The picture 

Let us begin with some definitions. 
A limit cycle of a renormalization group {TZs | s G 
with fixed a and is a periodic orbit of the flow Uic 



{[//;^ |s e K} where U^ 



(s+A) _ jj(s) 



Ic 



C/jj. with fundamental pe- 
riodicity A > 0. With a slight abuse of notation, we use 

(s) 

Uic to denote the orbit itself, while U^^ denotes a point 
on it. Let us call a function u discretely self-similar with 
parameters (a,/3) if it satisfies u{t,r) = u'^^'"'f^\t,r) 
for some A > 0. Each discretely self-similar solution 
u of the PDEs with parameters (a, 0) of the Assump- 
tion S in Sec. II B is related to a limit cycle Uic by 
u{t,r) = (-t)"C/^~*"(r(-t)-^). 

For a limit cycle C/ic, it is convenient to con sider eigen- 
modes of T. 



mode F of T, 



where Tg^jj is defined in (2.7). An eigen- 
(0) is a function satisfying 



for slightly supercritical solutions, where k is the eigen- 
value of the unique relevant mode of ^(o) . 

Here "self-similarity" in (SI) and (S2) should be un- 
derstood as "discrete self-similarity." 

Rough idea of the 'proof is as follows (see the following 
section for rigorous results). We reduce the problem of 
a limit cycle to that of a fixed pointf^ by considering 
the Poincare section S, a codimension-one submanifold 
transverse to the limit cycle, and the Poincare map V 
on it induced by the renormalization group flow TZ. (See 
Fig. |.) {S,V) deflnes a discrete dynamical system, for 
which the limit cycle in F corresponds to a fixed point in 
S. The behaviour of the map V on the Poincare section 
S is the same as that of the map TZa on the whole phase 
space F in the continuously self-similar case. Thus it 
follows from Claims 1 and 2 that Claims 1' and 2' hold. 



2. The Poincare map and the local behaviour of the flow 



■cA 



F. 



(2.94) 



The exponent k, which is equal to the Lyapunov expo- 
nent, gives the mean growth rate of the perturbation F. 
Relevant, irrelevant, and marginal mod es ar e defined ac- 
cording to the sign of Rck as in Sec. II A, It is easily 
seen that if F is an eigenmode of 7^ ^(o) then T jj(o)F 

is an eigenmode of 7^ jjU) with the same k. Note that 
for a limit cycle Uic there is a trivial marginal mode Fq 
of 7. rjio}F which is the tangent vector of C//°^ in F: 



_fn = lim 



u 



u 



(0) 



mic 



(2.95) 



To measure the 'distance' between a vector U and the 
limit cycle Uc, we introduce 



In this and next sections, we give rigorous results under 
certain assumptions. In this section, we consider the local 
structure of the phase space and show how the Poincare 
map reduces the problem of the limit cycle to that of the 
fixed point. 



^^By reducing in this way, we are losing some information 
on detailed behaviour of the flow. For example, if one is 
interested in a 'wiggle' predicted in [Q, one should study the 
continuous time evolution directly. 

^■^ After the submission of this paper we learned that the 
"wiggle" has been studied and numerically confirmed by Hod 
and Piran [Eel. 
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Let us first list assumptions used in the derivation, 
which are similar to those in the case of continuous self- 
similarity. To simplify the notation, we often write 7a 
for 7^ ^(0) in the following. 

Assumption LI' (Invariant subspaces of 7a) The tan- 
gent space T^(o) F of F at U^^'' is a direct sum of invariant 
subspaces of 7a: 



r^(o,r = i?"(f/f ) ® E'iU;-") e E%U,'^'); (2.98) 



(0)^ 



where E'^iul^^) is a one-dimensional eigenspace spanned 
by Fo satisfying 



(2.99) 



and F"(C//J'^) and -E"(J7j[!'') are of dimension N and of 
codimension A'^ -I- 1, respectively, and there exists k > 
such that 



||rA(F)|| >e«^||F||, FGF"(C/f ), 
Ta(F)|| <e-«^||F||, FeF^(C/if). 



(2.100) 



Assumption L2' (Smoothness of TZ.) The renor- 
malization group transformation TZ considered as a func- 
tion K X r 9 {s,U) ^ TZsU e r is of class C^'^ in a 
tubular neighborhood of the limit cycle U\c- Namely, for 
\s'~s\,\\U'-U\\ < (5o andsupo<,„<,dist(7^,-([/),C/lc) < 

Us'U - UM = {s' - s)n{nM) + Oils' - sp), 

(2.101) 

Tz,u' - = %,u{U' ^U) + 0{\\U' - Uf). 

(2.102) 

Remark. Assumption LI' together with the existence 
of TZ in Assumption L2' is a natural limit cycle version of 
Assumption LI of Sec. IIC. Assumption L2' is slightly 



stronger than Assumption L2. However, we consider it to 
be also physically reasonable and adopt it in order to rule 
out subtle and unexpected behaviour of the flow near the 
limit cycle U\c- Assumption L2' implies the following: 

TZ{U') -TZ{U) = 0{\\U' -U\\), (2.103) 
{%,^u - %,u) F = 0{\s' - s\ ■ \\F\\), (2.104) 
{%,u' - %,u) F - Oi\\U' - U\\ ■ \\F\\). (2.105) 

Let us define the Poincare map and reduce the prob- 
lem to that of a dynamical system on a submanifold of 
codimension one. 

Let S" be a submanifold of T defined by 



(0) 



Y, r = £:"([/, 



i?^(C/f ) 



(2.106) 



A point moving along a flow curve TZsU starting from 
point U on S will return to S for some s ~ A, and we can 
define the Poincare map T' : S ^ S so that P(C/) G S. 
More precisely, we have the following lemma. 



Lemma 3 (The Poincare map) 

Under Assumptions LI' and , the following (i) and 
(ii) hold. 

(i) There is a 5 -neighborhood of U^^'^ in S and the 
Poincare map T' which is a local C^'^ -diffeomorphism 
from S to S satisfying 



nU)^TZp^u)iU), 



(2.107) 



C^'^ and p{U^^^) = A. In particular, 



where p : S ^ is \^ uuu, pyji^ 
the tangent map T^qj of T' at ul^"^ is given by 



v: 



T, 



(2.108) 



(ii) Furthermore, a tubular 6 -neighbourhood of U\c is fo- 
liated by sections 



jY, (se 



(2.109) 



and on each section there exists the Poincare map 7^^*^ 
which satisfies all of (i) with Y, T' , S and ul^"^ replaced 
by Y^^\ T'^''\ S"*-*'-* and uj^K In particular, we have 

(2.110) 



The Poincare maps T''") for various s are related by 

V^"+''Kv,, =V,, oV^''\ (2.111) 

where Vs' ■ U ^ TZp^u){U) (s e I = [-A, A]; is 
a foliation-preserving local diffeomorphism on F which 
maps S'*-*-' to S^^~^^ K Here p is a C^'^ -function deter- 
mined by s' such that TZp(ij){U) G S^^^^ \ and satisfies 

p{u[^^) = s' . The tangent map T^f, ofVs o,t U\c ^ satis- 

fies 



T' 



Moreover T' is C^'^ , 
I X ((5— neibourhood of U\c)- 



(2.112) 

a function on 



Remark. The choice of S being a linear space, though 
in general it need not to be, makes the subsequent dis- 
cussion much simpler. In particular, by reducing the flow 
on F to maps on S, we can use the results in the fixed 
point case. 

Remark. In fact, (i) is included in (ii). One could define 
Vs first and then the Poincare map by 'P^*-' = 'Pa|s(=)- 

Proof, (i) The proof here is based on the standard treat- 
ment of the Poincare map (see, e.g. [|l], Section 8]). As- 
sumption LI' immediately implies that S is transverse to 
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?7ic at i.e., the direct sum of the tangent spaces is the where the inverse derivative of /, evaluated at 



whole tangent space, since T.-mS ^ £;"([//°^)®£;^(C//°^) {U,p{U),r{U)), exists and is given by ( |2.114D 

andT^(o)C/ic = £;=(t//,°^). 

Consider a map / :rxMx5^r, f{U, s, U') = 
TZsU - U'. We have f{ul°\A, ul°^) = 0. We show that 



and 



we can implicitly define maps p ; B{ul^\ 6q, T) - 

V : B{ul°\S'f;,T) Shy f(U,p{U),V{U)) = for some 

positive J^,', where B{u[^\5'i, T) is an open ball of center 

and radius 5q in F. 
Assumption L2' implies that / is C^'^, and the explicit 
form of the derivative of / at (f7, s, U') with respect to 
(s,t/') is 



The Poincare map on S is the restriction of to n 
B{u[^^ , 5q, r), which we denote by V. The transversality 
of TZ{U) an d S im plies that this restriction is a one-to-one 
map. Eq. (2.120) implies that 



di a/ 

ds ' dU' 



R X y 



tn{n,u) - F, 



(2.113) 



where df /dU' is understood as a Frechet derivative. If 
il{TZsU) ^ Y this has a smooth inverse 



a/ 5/ 

ds ' dU' 
F 



T^f{u,s,u')iF)\\ 



-^fiu,s,u')iF) (2.114) 



\n{nM)\\ 

for F G T/([/.s,(7')r, where ttu is the projection of T^/F 
into 5 along TZ CU). At (f//°\ A, [//°^) there exists a 



smooth inverse ( 2.114 ) because of the transversality of 
Uic and S, i.e., 7^(7^AC/lc) = ^^o ^ i"- 

Then the implicit function theorem implies that 
there exist Sq > and a one-to-one map {p, V) : 
B{u[°\5'iX) 3 U ^ {p{U),r{U)) gRx S satisfying 
f{U,p{U),V{U)) = 0, and that {p,V) is as smooth as 
/, i.e., of class C^'^. Denoting its partial derivative by 
{qu,T^), we have 

piU') - piU) = qu{U' -U) + 0{\\U' - Uf ), (2.115) 

V{U') - V{U) = T^{U' -U) + Oi\\U' - Uf), (2.116) 
iqu' ~qu)F^Oi\\U' ^U\\-\\F\\), (2.117) 

F = 0(||C/' - ;7|| • ||F||). (2.118) 



The explicit form of the derivatives are given by 



Qui-) 
which follows from 



mu)\\ 

u ° 'Fp{u),Uy 



dl df_Y 
ds ' dU' I 



(2.119) 
(2.120) 

(2.121) 



T, 



(2.122) 



where the second equality follows from invariance of Y 
under T^ jjm ■ 



r(0) 



(ii) It follows from Assumption 1' and T. (s) o T ^, 

j(o) E{ul^''ys are invariant sub- 



(0) that T j^.. 



spaces of 7^ jj(s} so that S^^^ and ii-iU^^ ) are transverse. 

Then by the same argument as (i), one can define Vs 
and the lemma is proved. (The existence of a tubular 
neighborhood on which Vs is well defined follows from 
compactness of /.) □ 

Though we will consider qualitative behaviour mainly 
on the Poincare section, we show here the qualitative 
behaviour in the whole phase space, which is the limit 
cycle version of Proposition 0(i). 

Proposition |l|' (Invariant manifolds of TZ at U\c) 
Under Assumptions LI' and LS! , we have the follow- 
ing: In a tubular neighbourhood ofU\c, there exist a sta- 
ble manifold W^{U\c) of codimension N, and an unsta- 
ble manifold W^{Uic) of dimension N -\- 1, whose tan- 
gent spaces at uj"^ are T ^(o) (-B"(f//°^) ® -B'=(f//°^)) and 

T^^(o)(£;"(t//°') e (E^ulf)), respectively. 

In terms of the Poincare sections, there exist a sta- 
ble manifold W''{ul'\s'-'^) and an unstable manifold 
iy"(C//;\S'(")) of each {S'^'^V'^''^) and are given by 



ly " (c//;^ , ) = Vs w''{ul°'> , s) . 



(2.123) 
(2.124) 



The stable and unstable manifolds W^(Uic) and W^(JJic) 
of the whole flow are written as 



w%Ui,) = [nsW%ul^\s) 
^^"(c/l,)-{7^,M/"({//^"^5) 



s G ] 

s e 



(2.125) 
(2.126) 



Proof. We now show that discrete dynamical system 
{S,V) satisfies Assumptions LI and L2 of Sec. II C, i.e., 
we show that the assumptions hold under the replace- 
ment 



(F,7^,,r,,^/)— ^(5,p,rJ'). 



(2.127) 
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Assumpt ion L2 holds f or (S', "P) because V is of class 
C^^i (see ( |2.115| )- (^.118|) ). Assumption LI holds for 
{S,V) because of and Assumption LI'. Then 

Proposition |l| of Sec. [I C implies that in 5* there ex- 
ists a stable manifold of codimcnsion iV in 5' and an 
A'^-dimensional unstable manifold of ul^\ Again, more 
sophisticated argument guarantees their smoothness p^ , 
Theorem 5.2, Theorem 5.n.4]; they are of class C^'^. 

Next, consider a one-parameter family of discrete dy- 
namical systems (S''^*\ P'''^). Lemma raii) states that 
Vs is a diffcomorphism from S to S^'^^ which pre- 
serves the action of the Poincare maps on them. This 



implies that lim„_oo = 



lim„_oo^"(P-.C/) = 
manifold W^iuf^ 
that 



(0) 



S) to W{ul''^ 



so that Vs 
Sis) 



if and only if 
sends the stable 
. It also implies 



is tangent to ^^{s)E^{U^^ ) so that 



,(.,(i?^(c//;')ei?^(f//;))). 



We 



WiU^^') is tangent to T 

have a similar result for the unstable manifold (consider 
'P~^ in place of V). Thus there exist the stable and un- 
stable manifolds W^{Uic) and W^{Uic) of codi mcnsio n N 
and dim ension iV-|- 1 in F, which are given by (2.125) and 
(2.126), respectively. □ 



We now consider quantitative behaviour of the renor- 
malization group flow. We assume the following, which, 
because of Proposition 0', asserts that the stable mani- 
fold is of codimcnsion 1 and the unstable manifold is of 
dimension 2. 

Assumption Ll'A (Uniqueness of the relevant 
mode) £;"(f//°^) is one-dimensional, i.e., = 1. Thus 
^"(^i"'') eigenspace spanned by F"^"' satisfying 



where k is the unique positive eigenvalue of 7a 



(2.128) 



Consider the time evolution U^^^^ = TZs{Uu 



initj of an 

initial data Uinit in Si''°\ We are going to trace the 
evolution on the Poincare section S. To measure how 
[/init is deviated from the critical surface, we introduce 
a decomposition of U-mit in the tubular neighbourhood 
of Uic- Lemma |^(ii) implies that for any U-mit there is 
a uni que can decompose any Uinit as 

(ra on 9 Uinif- 



{so) 



A^o) 



(2.129) 



As in the continuously self-similar case, we generically 
have 6(0) = 0{p — p^)- In the special case of t/init € S", 
Proposition § applies for the dynamical system (5*, 'P) 
and we immediately have the estimates in Proposition 0. 
In a general case where U-mit ^ S (say, the point A of 
Fig. 0), all we have to do is to estimate how U±_ and b 



behave before the flow intersects S for the flrst time (at 
the point B of Fig. Q); after that, we can directly apply 
Proposition ^. The estimate is given by the following 
lemma. 

Lemma 4 Under Assumptions LI' , LS! and LI' A there 
exist positive constants 5 and e which satisfy the follow- 



ing: Let \b\ < e, \\U± ~ Ul^°''\\ < S, s' 



So and 



U^U^ + bF['°\ [/^ e M^"(C/ic)nS'("''\ (2.130) 
U ^Vs'U ^U^ + bFi, U^£W''iUic)nS. (2.131) 



Then 



U±-Ul^^\\<C{S + e'), 



C 



<b<Cb, 



(2.132) 
(2.133) 



where C is a positive constant. 



Proof Let b' = b/WF^""^]] and b' = b/\\Fi\\. We have, for 
fixed So 7 

U ^Vs'iU^+b'F^''''^) 

= Ps'iU^) + b'T^l'F^''^ + 0{b'^) 



^VAU^)^h'T'i, Fr 



o 



{\b'\ {\b'\ 



(so) 



^PAU±) + b'F,+o{\b'\{\b'\ 



(so) 



(2.134) 



where we have used ( 2.112 ) at the fourth equality. {0{x) 
above may depend on s' or sq.) Since Vs>{U±) G 
Wiul^^), we have 



tj^^Vs'iUA + Oieie + S)), 
b' = b'{l + 0{e + d)). 



(2.135) 
(2.136) 



The first equality gives 

\\u^ - c/f II < WVAUA ~ vAv\:'H\ + o w + i)) 



0{6 



(2.137) 



where we have used the smoothness of V in Lemma ||. 
Eq. ( |2.136D gives 



b^b-\\F^'"^\\il + 0ie + 6)) = b-0il). 



(2.138) 



Since < s' < A, the 0{x) terms above are bounded by 
a constant independent of sq. □ 
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3. Global behaviour of the flow 

Now that the local flow structures have been analyzed, 
we can state a rigorous theorem which corresponds to 
Theorem |l|, under moderate assumptions on the global 
structures of the flow just as in the case of continuous 
self-similarity. One possibility is: 

Theorem 2 Suppose Assumptions LI' , LS! , LI' A, G2, 
and G3 hold, and Assumption Gl holds for the discrete 
dynamical system {S,V). Then the following hold. 

(i) All the initial data sufficiently close to Uic (but not ex- 
actly on W^{Uic) once approach, but finally deviate from, 
Uic- The final fate of these data are either a flat space- 
time or a black hole, depending on which side ofW^{Uic) 
lies the initial data. 

(ii) Consider the time evolution of the initial data ( 2.12(\ ) 

with \\U± - U['°\\ <5,Q <\b\< 5, which is on the side 
of black hole of W^{U\c). Then as b goes to zero, the 
mass of the black hole formed satisfies 



A/, 



c^\bf/''<^<c,\bf/- 



M, 



(2.139) 



with b-independent positive constants Ci , C2 . Here Mq is 
the initial (say, ADM or Bondi-Sachs) mass. 

Proof. Lemma ^ reduces the problem into the dynamical 

rstem (5, V) on the Poincare section S, where Theorem 
can be applied. □ 



E. Universality class 



problem reduce to the cases treated in previous sections. 
As a result, we understand why and to what extent con- 
dition (U3a) holds. Knowing the universality class (the 
class of models which exhibit the same critical behaviour) 
is very important from a physical point of view, because 
it guarantees that certain models exhibit the same critical 
behaviour with a model in the same class, whose anal- 
ysis could be much simpler. (Imagine we want to deal 
with a realistic matter. It would be almost impossible 
to know or deal with the exact form of the microscopic 
interaction. Universality could make the problem easier 
and accessible, by reducing it to that for a simpler model, 
e.g. perfect fluid.) 

Suppose unknowns u satisfy the PDEs of the following 
form, which is not necessarily scale invariant: 



, du du 



= 0. 



(2.140) 



Let us define the scaling transformation Sis, a, (3) : 
u{t,r) ^ «(")(*, r) = e°"'u{e-H,e-'^''r) as in {^). Then 
u^'^^ satisfies equations 

L^'Uu^'\^,^,t,r]=0, (2.141) 



dt ' dr 



where 



e'^'Li e-^'u^'i, e 



dt 



dr 



In this section, we focus on the problem of the uni- 
versality in the original sense of the term in statistical 
mechanics, based on renormalization group philosophy. 
We have largely benefitted from the formulation of [[l8[ . 
This is the rigorous version of the heuristic argument of 
universality class given in (l6| . |^ 

Consider general cases in which the PDEs do not nec- 
essarily satisfy the fundamental ass umption of scale in- 
variance. Assumption S in Sec. II B, or self-similar solu- 
tions with the scale invariance is not relevant for a criti- 
cal behaviour (e.g. having more than one relevant mode). 
Methods of renormalization group gives us a clear under- 
standing on which systems can exhibit the same critical 
behaviour. The main point is that the renormalization 
group transformation drives the equations of motion to a 
fixed point, where they gain the scale invariance, and the 



H,e 



(2.142) 



After the first draft of this paper was submitted for pub- 
lication, we have learned that a heuristic argument of univer- 
sality class has also been presented by Gundlach and Martin- 
Garcia , which is at the same level as had been presented 
in 0. 



and 7 = (71,72, ■•■,7m) (where m is the number of com- 
ponent of L) is a set of constants which may depend on 
a and (3 and are chosen so that i^*' remains finite and 
nonzero when s ^ 00. It should be noted that if L^"' = L 
holds with some 7 the system L = is in varia nt under 



scaling transformation in the sense of Sec. II A 



A renormalization group transformation is now defined 
to be a pair of a transformation on the EOM: L v-^ L'-"', 
and the following transformation on the F: 

7^s.L : U,{i) = u,{^l,0 

="!''(-!, e) = e"'%.(-e-^e-^^C), (2^143) 

where we added subscript L to show explicitly the RGT's 
dependence on L. See Fig. |[ We emphasize that the 
renormalization group transformation is now to be con- 
sidered as a pair of two transformations on different 
spaces: the space of EOM, and the space of initial data 
(phase space F). One can easily show the (semi)group 
property 



(2.144) 
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(a) 



(b) 



\L 




FIG. 5. Conceptual diagram of RGT: (a) the case V — Q, 
(b) the case with 'irrelevant' V . To simplify the figure, only 
the unstable manifold is shown. 



goes to zero. In fact, if one starts from the same ini- 
tial condition this usually does not hold (e.g. the critical 
value Pc{V) is shifted due to V = L — L*)\ we have to 
prove structural stability for these systems. 

In the following, we provide a sufficient condition un- 
der which the above naive guess does hold. To simplify 
the notation, we define F*^*^ = L^"-* — L* , and introduce 
1 1 |V^(C/) 1 1 1 which measures the "strength" of the difference 
V{U). The strength ||| • ||| should be so chosen that the 
following Assumption 2C should hold. We write the time 
evolution operator in terms of L by TiY , and that by L* 
by TZ] and their tangent maps by the same superscript 
convention. We assume that the finite time evolution op- 
erators are regular with respect to this difference. Con- 
cretely, we assume ( in ad dition to Assumption LI of and 
Assumption L2 Sec. II C, which are about time evolutions 
without V) : 

Assumption L2V (Regularity in the presence of 

V) There are positive Ki,K3,do and a such that the 
following hold for ||?7-C/*|1, < So and |||y(?7)||| < Sq: 



Kiu + F) = nZ{u) + r^jj{F) + o{k^\\f\\ 

\\nZ{u)-n„{u)\\<m\v{u)\\\ 

l!r;;^(^^)-r..t;(^^)||<X3llim)lll-ll^^ll 



(2.147) 
(2.148) 
(2.149) 



s = s„ > s„_i > ... > So = 0. (2.145) 

This states that long time evolution problem of the orig- 
inal equation defined by L is the composition of finite 
time evolution problems of equations with (in general) 
different PDEs = 0. 

We now turn to the problem of determining the uni- 
versality class. Concretely, we consider the following sit- 
uation. Suppose, with a certain choice of a, and 7, 
L*-"^ approaches a fixed point (i.e. a fixed function of its 
arguments) L* asymptotically: 



Lis) 



(s->oo), (L*) 



is) 



L* 



(2.146) 



It follows from definition ( p.l42| ) of transformation of L 
that the fixed point L* = satisfies the scale invariance 
in the sense of Sec. |ll A| . Suppose, in addition, there ex- 
ists a fixed point (or a limit cycle) of TZs,l* which has a 
unique relevant mode. And the question is: under what 
circumstances do L and L* exhibit the same critical be- 
haviour (i.e. in the same universality class)? 

At first glance the question looks rather trivial; a naive 
guess would be that they are always in the same univer- 
sality class, because L*^''^ approaches L*. However, there 
is a rather subtle problem involved here: Values of i^^^ 
and L* at the same argument will approach each other 
as s ^ 00; but this does not necessarily mean that their 
difference at their respective (generally different) [/'■''', 
which are defined by the time evolution by L^^^ and L*, 



These guarantee the structural stability of the flow in the 
neighbourhood of U* and V = 0. Concretely, we have 
the following Pr oposi tion, which corresponds to Proposi- 
tion! (i) °f S^c. H^. We denote by E%U*) and £;"([/*) 
the invariant subspaces of 7^,1/* , in the absence of V . 

Proposition 3 Suppose there is a fixed point U* of 
TZs,L' with a unique relevant mode, i.e. Assumptions LI, 
L2, LI A of Sec. hold, and \\\V^^\U)\\\ for 

\\U — U*\\ < Sq. Suppose further that Assumption L2V 
holds. Then there exist small positive constants S'^ , (5" 
(how small is specified in the proof) such that if 

ll|y(^)(;7)||| <(5;' for \\u-u*\\<5[ (2.150) 

then for any G E^{U*) with \\F^\\ < 6[, there exists 
F^"(i^") e £;"([/*) such that 



TZnaiU* +F^ + F^{F''))^U* 
s n — > 00. 
We define 

W''{U*;V) = {U* + + F^iF"")} 



(2.151) 



(2.152) 



as we did in (2.52). Note that now the definition depends 
on V. 

Sketch of the Proof of Proposition |^. The proof is done 
by Bleher-Sinai argument as was done for Proposition 
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We sketch the proof by mainly pointing out necessary 
modifications. 

Step 1. Decomposition. Same as the proo f of P ropo- 
sition |l|. We employ the same decomposition ( 2.37 ) , and 
for each F^^ try to find Oc such that 



(2.153) 



Step 2. Recursion. By (2.148), we have 

nX{u* + F) = n^iu* + F) + oiKsiwviu* + f)\\\). 

(2.154) 



Calculating ^^(;7*+F) .just as in ( |2.42| ), we this time 
obtain [cf. (|l43| ), (|]5|), ([is^ )] 



(2.155) 



F:,+,=%,U'{F^)+0{Ki\\Fr,r) + 0{K3Vn), (2.156) 
where u„ = sup 1 1 (t/)| 1 1. Note that the only 

l|C/-C/-||<<5o 



difference from Sec. [I C 1 is the presence of 0{K3Vn). 

Step 3. Continuity. Same as the proof of Proposi- 
tion |l|. Assumption 2C guarantees the continuity in ao, 
even in the presence of V. 

Step 4- Solving the recursion. Same as the proof of 
Proposition |l]. Our goal is to find sequences rn,tn for 
which (a) — (c) of Step 4 of the proof of Proposition |l| 
hold. 

Step 5. Proof of (a) — (c). This part significantly dif- 
fers from the proof of Proposition ^ , due to the presence 
of the term 0{K^Vn) in the recursion. This time, we de- 
fine tn = rn as follows. First define u„ = supj,>„ Vk and 
choose 5 1 and 5'{ as 



5[< 



l-e- 



1-e- 



-5[. 



(2.157) 



27^1 ' ' 2K3 
Then define r„ = i„ recursively, starting from tq = S'l as 
rn+i = (e-^" + KiS[) r^a + K^v^ (2.158) 



or (A 



KiS[ < 1) 



n-l 



(a) S[+K,Y,yk{A) 



2-fc-l 



(2.159) 



k=0 



We now show that (a) — (c) is satisfi ed by this choice 
of r„. First, dividing the sum over k in (2.159) according 
to , (2.159) implies {v^ is non-increasing in k by 

definition) 



1/2 



k=0 

n-l 

i:(A') 



/ \ 7i — k — l 



n — k — l 



(2.160) 



Three terms on the RHS goes to zero (as n oo), either 
because e~'"^ -|- KiS[ < 1 or Vn/2 0. This proves 
(a). Similar reasoning also shows the following uniform 
bound: 



< 



(A)"^l + i.3E^i'(A') 



/\ n — k—l 



< 2S[. (2.161) 



Now we turn to (b) a nd (c). First note that it is suffi- 
cient to prove [cf. (|2.60|), (|2.6lD] 



Tn+i > e ''"rn + Ki{r„ V t^f + K-iV„ 
tn+i < e'^'^tn - Ki{rn V tnf - K^Vn- 



Because we are taking r„ = tn and because r„ < 25'i, it 
is then sufficient to prove 



Tn+l < {e''" ~Ki5[)rn-K3Vn. 



(2.162) 
(2.163) 



Eg. ( 2. 162) is satisfied, by our definition of r„. And 
( [2.163 ) is satisfied, if 



2Ki5[)rn > 2K3Vn 



(2.164) 



It is not difficult to see (by induction) that (2.164) is in 



fact satisfied under our choice of 6[ and 6'{. This com- 
pletes Step 5. 



□ 



To complete the story, we have to make additional as- 
sumptions on the global behaviour of the flow in the pres- 
ence of V. Concretely we assume: 

Assumption GV (Global behaviour in the pres- 
ence of V) Consider the time evolution of initial data 
U± e ^±{€2,63) in the presence of V. Then there exist 
£2,53, 64 > such that if 1 1 |T^(J7±)| 1 1 < S4 then Assump- 
tions Gl, G2, G3 hold for suitably chosen C, C . 



In the above, A±{e2,S3) was defined in ( 2.75 ), in terms 
of dynamics without V. 

Under these preparations, we can now state our theo- 
rem. 

Theorem 3 Suppose there exist a, 7 and /3 > satisfy- 
ing L^''^ — > L* {s ^ 00), where L* is a fixed point of the 
scaling transformation: (L*)^'*) = L* . Suppose there is a 
fixed point or limit cycle U* ofTZs,L-' with a unique rel- 
evant mode, i.e. Assumptions LI, L2, LlA of Sec. 11 C 
hold. Moreover, suppose Assumption 2C and Assump- 
tion GV hold. Then, 

(i) if there exists > such that for \\U ~ U*\\ < Sq 
we have (1) \\\V^'KU)\\\ ^ as s ~> 00 and (2) 
|||y('*)(t/)||| < 5^ for sufficiently small 65, then the sys- 
tem L — exhibits qualitatively the same critical be- 
haviour as the system L* = 0. (i.e. there is a critical 
behaviour and near-critical solutions once approach and 
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then deviate from U* ). 

(a) If in addition ~^ sufficiently fast (for 

example, integrable in s on (l,oo) is suffi- 

cient), the system L — exhibits quantitatively the same 
critical behaviour as the system L* = (i.e. the same 
critical exponent). 



Remark. In Assumption GV, V in fact means 
at the instant when the flow exits the perturbative re- 
gion through the tube A±{e2,S3): the local analysis 
which corresponds to that of Proposition ^ shows that 
the near (but not exactly) critical solutions once ap- 
proach and then deviate from U*, and they exit from 
the perturbative region through the tubes A±. Because 
|||P"(")(C/)||| ^ as s ^ oo, the effect of V goes to 
zero as p ^ Pc{V), and could be neglected. However, 
we made the Assumption GV in order to avoid patho- 
logical situation where the slight effect caused by this 
small leads to qualitatively different global 

behaviour. (Because the solution is expected to blow up, 
we cannot expect nice regularity in V in this region.) 

Theorem ^ explains how and why claim (U3a) of 
Sec. II B holds. The mere condition L ^ L* is not suffi- 



cient to guarantee the universality. The following proof 
in fact suggests the deviation in the critical behaviour: 



Mbh(x b-pc|^/"/(p) 



where 



f{p) « exp 



C2 



;i log \p-pc 



V{s)ds 



(2.165) 



(2.166) 



with positive constants ci,C2. Concretely, f{p) will be a 
power of logarithm in |p — Pc| for IHV^**)!!! oc 1/s, and 



f{p) cx exp (const, log^ " b ~ Pc|) 



(2.167) 



for |||F(^)||| cx (0 < a < 1). The deviation in { ^.Wt l 
is greater than logarithmic, although it does not change 
the critical exponent (if we define it as the limit of the 
logarithmic ratio of both hand sides). 

Proof of Theorem |^ Proposition ^ shows the existence 
of a critical solution. We now trace the time evolution of 

off-critical solutions. 

We employ the scheme of Sec. [ICS. That is, we de- 



compose the initial data U as [cf. (2.77)] 
and for n > decompose [cf. ( ^.78| )] 

7^L(f/) = 7^L(c/c) + J^. 



and further decompose F„ as in ( 2.79 ) 



(2.168) 



(2.169) 



(2.170) 



Note that Uc here means the critical solution in the pres- 
ence of V (and its ti me evolution). Then ar guing as we 
did in deriving ( 2.85 ), we now obtain, using ( |2.147 ) 



7^|;+l),(c/)=7^r(t/("'^) + i^„) 

V 



= 7^r(t/i" 

which implies 



)+T'^,..AFn) + OiK,\\F^r) (2.171) 



Fn+l — T 



,{F^)+om\Fr). 



(2.172) 



Now using ( ^.149| ) and ( ^^ 



T^^uiF) = [r^,u - (F) + [%,u - %,U'] {F) 
+ %.U'{F) 
= %,U'{F) + 0{K,\\\V{U)\\\-\\F\\) 

+ 0{&Ki-\\U -U*\\-\\F\\) (2.173) 



So decomposing F„ into F^ and F^ in ( 2.172 ), we have 
recursions for 6„ and /„ = 1 1 -F"^ 1 1 , gn 



(na) 



bn 



+ 1 



K^{\bn\y fn)0{Vn) 



U* 



(2.174) 



+Ki{\b,,\ V fn)[0{\bn\ V /„) + 60((7„)] 
fn+l < e '"'^ fn 

H\K\ V /„) [i^i(|6„| V /„) + 6i^i5„ + i^3«n] (2.175) 



As seen, the only difference here from (B3) and (B4) 
is the appearance of w„. The recursion is of the form 
considered in Lemma ^ (6i^if/„ -|-if3W„ playing the role of 

there), and this proves the quantitative estimate which 
corresponds to Proposition ^ In particular. Lemma ^ 
guarantees that the presence of 0{vk) does not affect the 
critical behaviour. 

The above concludes the local analysis of the flow in 
the neighb ourh ood o f [/*, which corresponds to that in 
Sec. [I C 1 and II C 3 The analysis of global behaviour, 
corresponding to that of Sec. II C 4, can be carried out in 



a similar manner under the Assumption GV. 



□ 



Remark. By refining the the proof of Proposition ||, we 
can show in generaf^ (Pc(0) denotes the critical value 
without V) 



^* Heuristically, we can go further. If V is of definite sign, 
we can heuristically get a lower bound on the difference, and 
show 



\p.{v)~pM\ = o{\\\v{u*)\\\). 



(2.176) 



This estimate gives an answer to the following question con- 
sidered in [Q: fix the value of p to the critical value when 
V = Q (i.e. p = Pc(0)), and vary the magnitude of V . What 
kind of critical behaviour can we observe? The above esti- 
mate, together with Theorem ^ shows that we will observe 
the same critical exponent /3bh = because 
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\Pc{v)-pM\<o{\\\v{un\\\). 



(2.178) 



We conclude this section by giving two examples to 
which the above analysis can be applied. Note that the 
following arguments are not rigorous in contrast to what 
we have been doing, in the sense that no rigorous verifi- 
cation of Assumptions are given. 



1. Example: perfect fluid with a modified equation of state 



of the fixed point, by simply taking |||V||| = /(p). Then 
Theorem || implies that the model belongs to the same 
universality class as that of / = 0, as long as (for fixed 

x) 



oo 



(2.184) 



This is the sufhcient condition we are after. An example 
of / which shows the same critical behaviour as / = is 



given by f{x) 



with < S < 1. 



The first example is a perfect fluid, with equation of 
state (/ is a given function) 



p= (7-l)p + /(p). 



(2.179) 



The EOM for this system is given by ( |3.3| ) of the next 
section, supplemented by the above equation of state. 
When / = 0, this is a usual perfect fluid, to be studied in 
detail beginning with the next section. It will be shown 
that the system with / = exhibits a critical behaviour 
governed by a self-similar solution U*, which is a fixed 
point of a renormalization group transformation induced 
from a scaling transformation S whose effect on p and p 
is: 



P 



'^l^'p{e~H,e-^'r), 



(2.180) 
(2.181) 



where /3 > is arbitrary due to gauge invariance. The 
EOM ( |3.3[ ) is invariant under the scale transformation, 
and so is the equation of state with / = 0: Lg = p — (7 — 
l)p. (Here we write Lg because this is the sixth of EOM.) 
Thus the syste m wit h / = i s an example to which our 
scenario in Sec. II B and Sec. QC applies. 

Our question is to determine the models (forms of /) 
which belong to the same universality class as the / = 
model. Now under the scaling transformation the EOM 
in question, Le = p— (7— l)/3 — /(p), is transformed into 



(^„i)p(^)_/(^)(p(^)) 



where 



fie 



(2.182) 



(2.183) 



At the fixed point in question, p* and p* = (7 — l)p* are 
finite smooth functions of their argument (Sec. IV). So 
we can expect Assumption 20 to hold in a neighbourhood 



Mbh{pc{0);V) (X \p,{0) ~ p,{V)f^^ oc |||V(I7*)|||''=^ 



(2.177) 



where we used Theorem ^ in the first step, and (2.176) in the 
second. 



2. Example: scalar field collapse 

Let us present another example, i.e. spherically sym- 
metric collapse of a scalar field. 

We consider a scalar field with energy-momentum 
tensor 

Tab = ^aq^^b^ - ^gab5'''Ve0V<i(/. - 5afcV^(0), (2.185) 

where gab is the metric tensor, V is any, say, C^-, func- 
tion, and the abstract index notation is used. The line el- 
ement of a general spherically symmetric space-time can 
be written as 

ds^ = —g{u, r) g{u, r)dv?' — 2g{u, r)dudr + r^dil^ . 

(2.186) 

We fix the gauge by the condition g{-, 0) = 1. The equa- 
tions of motion can be written in terms of retarded time 
u, area radius r and field variable h{u,r) as = 0, 

where functional L labeled by V is given by [|| 

Lv[h] ^Dvh^\{h- {h)) + \rg[h] V'{{h)), 

(2.187) 

where (•), 5, g, and Dv are functionals defined as 



drf, 



5[/]^exp(-4^r(r-i(/-(/))2)), 
5y[/]^(5[/](l-8^r2y((/)))), 

Dvf 



df i_ ..a/ 



(2.188) 



In the equation of motion, = above, met- 

ric variables g and g have been solved as function- 
als g and gy, respectively, of h, and the scalar field 
is given hy (f> = (h). Hereafter, for a functional F, 
F[h] denotes a function determined by a function h and 
F[h]{a,b) denotes the value of function F[h] at (a, 5). 
For functionals F — {■),g,gv,Dv, we rewrite the values 
F[h] (e^^u, e^^'^r) in terms of ft,'"-' (w, r) using the relation 
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{h){e-'u,e-'^'r) = e-"'{h^'^){u,r), (2.189) 

— 2as 

g[h]ie-^u,e-^^r)^[g[h(^^]y {u,r), (2.190) 

gy[/i](e-^u,e-^V) = ((g[M^)]) 

l-87rr2y('')(e-"''(/i("')))\ (w,r), (2.191) 



(2.192) 



\ ou 
dr 



(u, r) 



where 



y(-)(a) = e-2^"y(a). 



(2.193) 



We find the renormahzed functional (L y)'^^^ wi th pa- 
rameters a, (3 and 7 by substituting ( |2.189| )-( ^.192 ) 
into e'^^LY[h]{e~''u,e~^'^r) and reinterpreting it as 
(Ly)('')[/i(")](u,r), the value of function (Lv)(") at 
{u,r). 

Let us find fixed points, i. e., Ly satisfying (Ly)*-*' — 
Ly- It follows from ( 2.192|) -( |2.192D that we mu st choose 
a = 0, /3 = 1. Then (|2.189| )-( RT92|) and (|2.193D read 

(2.194) 
(2.195) 
(2.196) 
(2.197) 

(2.198) 



g[h]{e~''u, e'^r) 
gv[h]{e-''u,e-''r) 
{Dyh){e~'u,e''^'r) 
and 

Transformation of L is thus given by 

(Ly)(^)[/i(")](u,r) = e^^-''>'L[h]{e-'u,e-'^'r) 
= e(^-i)^Ly(.,[/i(^-)](u,r). 
Choosi ng 7 = 1, we have (Ly)'"-' = L,, 



5[M^)](u,r), 

5y(„[M^)](u,r), 

e-^Dyu,h^'^){u,r) 



-2s 



y(a). 



(2.199) 
, and from 

( ^.198 ) we find that Lq is the only fixed point, i.e., Lq — 
is the only model having scale invariance, among the 
models considered here. 

Choptuik's original result which was further con- 
firmed by a theoretical work by Gundlach IHJl^, was 
that the system Lq ^ shows critical behaviour with 
this choice of parameters, a = 0, /3 = 1, for the RGT. 
The system has a limit cycle which should have a unique 
relevant mode, according to our scenario in previous sec- 
tions. 



It follows directly from (2.198|) and (2.199) that 



^V* =0, (Ly)(^) Lo (s -> 00), (2.200) 

and the convergence is so fast that Theorem ^ is appli- 
cable. This shows that potential V is an irrelevant term 
which vanishes asymptotically. Therefore, all Ly's are 
in the same universality class as Lq, i.e., all scalar field 
models with arbitrary smooth potential V exhibits the 
same critical behaviour with the same critical exponent 
as the minimally coupled massless scalar field. 



III. EQUATIONS OF MOTION. 

We now begin our concrete analysis of gravitational 
collapse of perfect fluid, as a case study of renormaliza- 
tion group ideas. 



A. Equations of Motion 

The line element of any spherically symmetric space- 
time is written as 

ds^ = ~a^{t,r)dt^ + a^{t,r)dr'^ +r'^{d9^ +siTi^ Odcj)'^). 

(3.1) 

We assume the matter content is a perfect fluid having 
energy-momentum tensor Tab = pUaUb + p{gab + "aWb), 
where p is the density, p is the pressure, and is a unit 
timelike vector whose components are given by: 



Ut 



-a 



aV 



VI - "i^^ ■ 



(3.2) 



Here V is the 3-velocity of fluid particles. Introducing 
locally U = V/^/l - and W = 1/Vl - the equa- 
tions of motion (EOM) are given in this coordinate by: 



a r 1 — a 



a 2r 



+ A7rra^{pW^+pU^) 



a 2r 



— — ~A'Kraa{p + p)UW 
a 



+ AT:ra%pU^ + pW^) 



(3.3a) 



(3.3b) 



(3.3c) 



(apiy),+p(aT4^),+ ('^'"^^)-+/(^'"^)'- ^0 (3.3d) 



{apU).t + p{aU)^ + {apW)^r + p(aVK) -r = (3.3e) 

In the above, first three are obtained by taking linear 
combinations of (0,0), (1,1), (0,1) components of Ein- 
stein equation. The last two are obtained by taking lin- 
ear combinations of the "Bianchi identity" T'"';^ = 0. 
We have to supplement the above by the equation of 
state for the fluid. We consider the case 



P={l- 1)P, 



(3.4) 



where 7 G (1, 2) is a constant (adiabatic index). 

In terms of variables s = — ln(— t),a; = ln(— r/t), and 
introducing 



N=—, A = a^, u! = Anr'^a'^p, 
ae^ 



(3.5) 
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we can write the equations of the system in an au- 
tonomous form, which makes the scale invariance of the 
system transparent: 



A- 



2lo (1 + (7 - 1)1/2) 



1 - 



^ = -2 + A-(2-7V 



3(2- 



A 

lVV,s 

1 - y2 

7) 



A,^ _ 



(3.6a) 
(3.6b) 
(3.6c) 



+ {1 + NV)— + 



NV- 



2 + 7 



ANV + (2 - -f )NVLj (3.6d) 



W i — V ^ 

= -(7-2)(7-l)7Vtj + ^^^^^A^+^-^AiV (3.6e) 

Only four out of the above five equations are inde- 
pendent. More precisely, the eqn. ( 3.6c ) is au to mati - 
cally satis fied by solutions of the set ( 3.6a ), ( 3.6b ), ( 3.6d ) 
and (3.6e), as long as they satisfy a boundary condition 
A{s, —oo) ~ l,y(s,— 00) = a;(s,— 00) = 0. In view of 
this, we pick up the above set of f our e quations as our 
basic equations of motion, and use (3.6c) as an auxiliary 
equation at appropriate stages in the following. 



IV. THE CRITICAL (SELF-SIMILAR) 
SOLUTION. 



To carry out the first step of our scenario, we first 
have to find out fixed points (i.e. self-similar solutions) of 
RGT. In this section, we find out (almost) all of the self- 
similar solutions of the EOM. Because this paper is in- 
tended to present the "renormalization group" approach, 
we here mildly try to exhaust all possible self-similar so- 
lutions. 



A. Equations for self-similar space-time 

We first require that the space-time is self-similar, i.e. 
that N and A depend only on a;: N — Nss{x),A = 
Ass{x). Conversely, it can be shown (see Appendix 
??) that any spherically symmetric self-similar space- 
times can be expressed in that form if a freedom of co- 
ordinate transformation ( |3.7[ ) is used. Then it follows 
from { \i.6\ ) that Wss and Vss are also functions of x only: 
uj = ujss{x)^ V = Vss{x). In this sense, the space time we 
are interested in is a fixed point of the RGT. 

The equations for a self-similar solution, TZ{U*) = 0, 
are then given by omitting terms containing derivatives 
of s in (3.6): 



A^x 



I- A 



2uj (l + (7-l)y2) 

1-^2 



^ = -2 + A-(2-7)c. 



(4.1a) 



(4.1b) 



B. Gauge degrees of freedom 

The only coordinate transformation which preserves 
the form of (3.1) is 



which corresponds to 

(s,x)k^ (ri(s),x-,s + /-i(s)), 



(3.7) 



(3.8) 



where F-^{t) = -e"-^ with / = dfjds ^ 0. Under 
this transformation, the variables h transform to h, where 



h{s, x) 



h{s,x) {h^A,u,V) 
f{s)h{s,x) (h^N) 



(3.9) 



Eqs.(3.6) are of course invariant under this transforma- 
tion. Fixing the value of at a point, for example, x = 0, 
of each constant s line determines the coordinate system 
completely. We shall retain the degree of freedom for the 
time being. 

One note: The gauge degree of freedom present in our 
syst em e nables us to take /3 = 1 in our formulation of 
Sec. [I A, without introducing new variables. 



(l + NV) 



A^^ _ 2'yNVu} 

~A ,x ~ ~ 1-1/2 



(4.1c) 



3(2-7) 



1- 1/2 

2- 1-7 



^AF -I- (2 - 7) Al/w (4.1d) 



b-l)(A'H-V,i^ + 2(l±^ 

= (2-7)(7- l)Atj + ^^^^--^A+ ^— ^AA (4.1e) 



from (4.1a) and (4.1c), one obtains an algebraic iden- 
tity 

1 I 2^(1 + (7 -1)^2) ^ 2^NVoj ^ ^^ 2) 



1-1/2 



1-1/2 



One could eliminate one variable (e.g. A) from (4.1) 
using the above identity. However, to avoid making 
equations too complicated, we keep using four variables 
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{A, N,llj,V) in the following, using the above (4.2) as a 
check at appropriate stages of our numerical calculation. 

U nder the self-similar ansatz, the coordinate freedom 
(3.7) reduces to 1 1-^ t — kt with constant fc, which corre- 
sponds to the translation of x [the other transformations 
alter the constant x lines in (t, r) space] . This freedom of 
coordinate transformation allows one to adjust the value 
of at a given point arbitrarily. We make use of this 
freedom, and fix the coordinate system by requiring that 
the sonic point (see below) be at a; = 0. 



A direct consequence of the Lemma is that any so- 
lution, which is regular at the center (i.e. A — > 1 as 
X — * — oo), must satisfy A> 1 for all x G M. 

Proof of Lemma ||. Follows from ( 4.1a ); The last term 
on the RHS is nonnegative, so we have 



A^> l~ A. 



(4.3) 

□ 



B. Conditions on self-similar solutions 

The behaviour of self-similar solutions has been exten- 
sively discusse d by Bogoyavlenskii in Q, followed by 
other works [pip^ ]. We here briefly outline our analysis, 
in our coordinate system, for completeness. 

We are interested in self-similar solutions which satisfy 
the following two properties: 

(i) The self-similar solution is analytic (or at least 
smooth, in the sense of having an asymptotic ex- 
pansion to all orders) for all a; g R, 



(ii) The space-time and the matter are regular, A 
and = 0, at the center {x = —oo). 



1 



The reason why we require these may be summarized 
as follows. In this paper, we are interested in self-similar 
solutions which could represent the space-time structure 
and radiation fluid profile at the edge of the formation of 
a black hole. Such a critical space time is formed from 
a suitable initial condition through time evolution given 
by the EOM, and we have introduced the self-similar co- 
ordinate (s,x) in order to absorb any singularity in the 
process. Thus, the solution should be smooth for all x. 
This is the first condition. 

The second (boundary) condition is a physical, natu- 
ral one; ^ = 1 at the center guarantees that the space 
time has no physical singularity at the center. V = 
at the center should hold for any spherically symmetric 
problem, except for the case where there is a source of 
fluid at t he ce nter. 

Now, (4.1) is a set of ODE's for four variables 
{A, N,LJ,V), which satisfies the Lipschitz condition ex- 
cept at the so-called sonic point (explained in sec- 
tion IV C). So, there are three types of possible singular- 
ities for (4T): a; = — oo, a; = oo and the sonic point. We 
consider each of them in the following, but before that, 
we conclude this subsection with the following simple ob- 
servation, which will be repeatedly used in subsequent 
sections. 



Lemma 5 Any solution of (^jj is monotonically non- 
decreasing, as long as A <1. 



C. Conditions at the sonic point 



1. The sonic point 



Eqs. (4. Id) and (4.1g) can be written in the form: 



a h 
c d 



(4.4) 



where a, 6, c, d, e, / are functions which do not depend 
on derivatives of ld and V . Solving the above in favor of 
Lij_x and V^x, the resulting equation violates the Lipschitz 
condition at the sonic pointy where the "determinant" of 
the coefficient matrix oi ui x and V x vanishes: 



det — ad — be 

^{[l + NVf-{^^l){N + Vf} 

« = 0- 



(4.5) 



Physically, the sonic point is characterized by the fact 
that the velocity of fiuid particles seen from the observer 
on the constant x line is equal to the speed of sound 



2. Expansion at the sonic point 

Because generic solutions at the sonic point can be 
singular, the analytic condition (i) required at the sonic 
point places a severe constraint on possible form of solu- 
tions and in the end leaves only one free parameter, as 
we now explain. 

The sonic point is characterized by the vanishing of the 
determinant, (4.5). In order to have finit e de rivatives, 



Lj^x and V^x^ at the sonic point, the rows of (4.4) must be 
proportional to each other: 



af — ec = 0. 



(4.6) 



These tw o co nditions, ( [4.5| ) and (4.6), together with the 
identity (4.2) between A, N,u!,V, enable us to express 
values of A, N, lo at the sonic point in terms of that of 
Vq = Vss (we write Ao,iVo,Vb, etc. for values of these 
functions at the sonic point). 

To find out the explicit forms of these expressions, we 
first list up the following properties to be satisfied at the 
sonic point: 
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1- -^0 > 0, because of its definition in terms of lapse 
function. 



2. Aq > 1, from Lemma |^ of Section [VB 



3. The derivatives uj^x and V,x should exist at the sonic 
point. In other words, first order expansion coeffi- 
cients, in particular those of lu and V must exist as 
real numbers. 



Note that ( |4.5| ), being quadratic in A^O; in general al- 
lows two solutions for A^o- However, one of them leads 
to < 1, and is excluded from the second requirement 
above. The first and third requirements together restrict 
allowed values of Vq to — 1 < Vq < ^7 ~ 1- -^-S a result, 
the zeroth order expansion coefficients are given by: 



No 



An 



LOq 



(4.7) 

72 + 47 - 4 + 8(7 - 1)3/2^/0 - (37 - 2) (2 - j)Vo^ 



72(1 -Fo2) 



72(1 -T/„2) 



(4.8) 
(4.9) 



Once these values at the sonic point are given, we can 
then determine the higher coefficients of power series ex- 
pansion of A,N,Lu, an d V with respect to x using the 
above conditions and (4.1). To avoid making presenta- 



tion too complicated, we do not give explicit expressions, 
but only remark that the equations for the first order ex- 
pansion coefficients, given in terms of Aq, No,luo, Vq, have 
two solutions in general; some of which are not allowed 
for some values of Vq due to the third requirement above. 

The self-similar solution which is regular at the sonic 
point is thus characterized by a single parameter Vq = 



D. Asymptotic behaviour of solutions as a; — > ±cxd 

In searching for desired self-similar solutions, it be- 
comes desirable to understand the correct asymptotic be- 
haviour of solutions as a; — !■ ±00. 



1. Behaviour at the center, x = ~oo 



Behaviour of solutions of ODE's (4.1), as x ^ —00, can 
be most easily understood by introducing a new variable 



M = NV 



(4.10) 



and rewriting (4.1) in terms of A, M,u;,V . It is easily 
seen that 



^= 1, 



M = , 

37' 



OJ ■ 



V ^0 



(4.11) 



is a stationary (fixed) point of ( [4.lD , i.e. all the x- 
derivatives vanish. Moreover, linear analysis around this 
stationary point shows that there is one unstable direc- 
tion for this stationary point, namely in the direction 
of AI + Thus, a generic solution which satisfies the 
analyticity condition at the sonic point will diverge as 
X —00, and only a small number of solutions with 
carefully chosen values of Vq can satisfy the regularity 
condition at x — ~oo] thus leaving only discrete accept- 
able solutions. 

By similar reasoning, it is easily seen that any self- 
similar solution which satisfies our boundary condition 
at the center (i.e. has a smooth derivative), allows a 
(possibly only asymptotic) power series expansion in ^ = 
at ^ = 0, and it is easily seen that the asymptotic 
behaviour is 



A^,^,{x) ~ 1 + A_ 



(4.12) 



with coefficients A_oo, -/V-00, w_oo, V-oo satisfying the re- 
lation: 



A. 



2lo. 



2 



(4.13) 



2. Behaviour at x = 00 



It is easily seen that 



TV = 0, A=l 



2tj{l-H (7- 1)^2} 

1-^2 



(4.14) 



is a stationar y po int (parameterized by two free parame- 
ters uj, V) of (4.1). Linear analysis around the stationary 
point shows that this stationary point becomes an at- 
tracting node for 



UJ < 



y2 



7(1 + ^^^2) 



(4.15) 



As will be explained in Section IV F, most of the self 



similar solutions which are analytic for a; < cxd found 
numerically satisfy this condition, and thus are seen to 
be attracted to this stationary point. 



E. Numerical Method 

Based on considerations presented in the above, we 
proceed as follows in order to find out desired self-similar 
solutions. 
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We first fix the value Vq at the sonic point a; = , 
and use the power series expansion there (Section |vcj) 
to derive the derivatives. Then sta rting from the sonic 
point X = 0, we solve ODE's (O) towards x = —oo, 
using Runge-Kutta fourth order integrator, with reliable 
error estimates. If at some x < (1) either A < 1 or 
(2) det = 0, we stop solving and conclude that this Vq 
does not give rise to a desired self-similar solution. The 
first ca se (A < 1) is excluded based on the Lemma || of 



section IV B . The second case implies the existence of 
another sonic point between a; = and a; = — oo. Such 
a solution with more than one (but with finite number 
of) sonic points can always be found by searching for its 
part between its left-most sonic point and x = — oo, and 
then solving (4T) from the left-most sonic point to the 
right (towards x = oo). However, it is quite unlikely 
that such solutions with more than one sonic points ex- 
ist, because the existence of only one fixed point already 
reduces the freedom of parameters to discrete ones, and 
other possible sonic point will further reduce it. 

Our numerical results show that each allowed Vq ends 
in either of these two possibilities. Moreover it is seen, 
both numerically and mathematically rigorously, that a 
desired self-similar solution exists at those values of Vq, 
where Vq -I- leads to the case (1) and Vq — leads to the 
case (2). 



F. (Numerical) results 

We have searched for self-similar solutions which sat- 
isfy the above conditions (i) and (ii) numerically, for all 
the allowed values of Vq. As has been extensively studied 
in there exists a sequence of self-similar solutions, 
each of which is characterized by the number of zeros of 
V. The number of zeros is an integer starting from 1, 
and the one with exactly one zero is the solution cited 
in 1^, which we call the Evans-Coleman solution. Due 
to a numerical difficulty, we could find only self-similar 
solutions with odd number of zeros of V. 

In Fig. 1^, we show the profile of the Evans-Coleman 
solution, which will play the central role in later sections. 
This is the same as Fig. 1 of |^ except for a which is due 
to the coordinate condition for t. 



V. PERTURBATION 

We studied perturbations around self-similar solutions 
by two different methods. The first one is to dir ectly find 
out eigenmodes by solving the system of ODE's ( 3.13| ) be- 
low, which characterizes the eigenmodes. The second one 
is to apply the so-called Lyapunov analysis. The former 
has the advantage that it reduces the problem to that of 
solving ODE's, and enables us fairly accurate estimate 
of eigenvalues, but has the disadvantage that it is almost 
impossible to search for every possible eigenvalues. The 




FIG. 6. Profile of a self-similar solution 
(the Evans-Coleman solution). Curves represent Ass (solid 
line), ln(A''ss) (dashed), ojss (dotted), and Vss (dot-dashed). 

latter has an advantage that it is quite suited to finding 
out eigenvalues in descending order of their real part, but 
has disadvantages that (1) eigenvalue estimates are not 
so accurate, and (2) we have to solve a PDE. However, 
we here emphasize that the PDE we have to solve for 
the Lyapunov analysis is a very regular one, and can be 
handled by standard techniques for solving PDE's. 

In this section, we explain the first method, that is, 
direct search for perturbations by shooting method. To 
simplify the notation, we express perturbations in terms 
of perturbations oi A = \n{A),N = \n{N),oj — ln(ct;). 
For example, the A-version of (p3h in fact means 



A{s,x) = Ass{x) -I- eAvar(s,a;) 



(5.1) 



where 



Ais, x) = \n{A{s, x)), A,,{x) = ln(A,,(x)). (5.2) 



A. Perturbation equations and the gauge degrees of 
freedom 

Perturbation eq uati ons are obtained by taking the first 
order variation in (^^) from the self-similar solution Hss 
considered: 



h{s,x) = Hss{x) + e/ivar(s,2:), 



(5.3) 



where h represents each of {A, N, cu, V). In the following, 
to simplify the notation, we write V for Vss and write 
Vvar for perturbations. 
Explicitly, we have: 



/O 



As 

\0 Cs Ds 



10 
10 
A^ 




dx 



I Ayax 
^var 
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/ Gi G3 G4 

Hi H3 

El E2 -B3 Ei 

\ El E2 E3 F4 



( ^var 
^var 

V Kar 



(5.4) 



with coefficients given in terms of self-similar solutions 
as 

7F 



7 



f-y2' 

7(iV + ^) 

1-^2 ' 



(5.5) 



C.^{^^\){N^V\ D.^ ^\^_yV ^ (5.6) 

(5.7a) 
+ (2 - i)u:NV 

(5.7b) 



Si 
^2 



(7 + 2) 



2 

6-37 
2 



NV 



ANY, 
2 + 7 



2 

1-1/2' 

= (2 - 7)cjA^T/, 
6 - 37 2 + 7 



Ea = 



-N- 



AN+{2- -f)LuN 



-NuJx 



-f{l + 2NV + V^)V,^ 
(l_T/2)2 



(5.7c) 



(5.7cl) 



El . '-^AN, 



(5.8a) 



7t - 6 2-37 
F2 = (2 - 7)(7 - 1)cj7V + -^^N + -—^AN 



1 -y2 ' 

i^3 = (2-7)(7-l)^A^, 

-/{N + 2V + NV^)V,^ 



-{j-l)Nu, 
2-7)(7- 

-F4 = -(7-lH 



Gi = 

G4 = 



-A, G3 

(1-^2)2 



(1-^2)2 

2{1 + (7- 1)F2}w 



(5.8b) 
(5.8c) 
(5.8d) 



1 -y2 



Hi = A, ff3 = (7 - 2)c^ 
We also have 

2-fNVLJ 



(5.9) 
(5.10) 



1-^2 



{Nvar + Wvar) 



27iVtj(l + F2) 



(1-V^2) 



2^2 



(5.11) 



We note that there is a gauge freedom due to t he co- 
ordinate transformation of order e, namely (3.8) with 
/(s) = s + e/i(s). The transformation ( |3.9| ) now becomes 
[taking 0(e) terms; ' = 



/ivar(s,a;) -f /i(s)/i^,(x) 

(^,<:^,"l^), 

'^var 

(s, x) + fi{s)h'g^{x) + fi{s) 

(5.12) 



This in particular means that one can always require 
^var(s, xq) = for a fixed xq. 



B. EOM for eigenmodes 

We consider eigenmodes of the form hvai{s,x) — 
hp(x)e'^^, with k S C being a constant. Substituting this 
form into (5.4) yields a set of linear, homogeneous first 
order ODE's for {Np, Ap,u!p,Vp). Exphcitly, we have 



/ 1 

10 

A^ B. 

VO G, D 




/O 



As B, 

Vo Gs D 



Gi G*3 G4 

Hi Q H3 

El E2 E3 E4 

El E2 F4 

A,. 
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(5.13) 



and 



kAp -)- dxAp — 



2-fNVuj 



r2\ 



2'^Nuj{l + 

(1-^2)2 



(5.14) 



Eq. ( 5.14 ) and the first row of ( 5.13| ) provide an algebraic 
identity : 



V 1 -1/2 

/M1+ (7-1)^^ 



l-y2 

/27tj{iV(l + l/2) + 2y} 

(1-^2)2 



lNV}\ _ 

1 Wp 



^P = 0, 



(5.15) 



just as in the case of self-similar solutions [cf. (4.2)]. As 
in our treatment of self-similar solutions, we here treat 
Ap, TVp, ujp^Vp as four unknown functions, and use ( 5.15 ) 
as a check at appropriate steps of computation. 

We require 7Vp(s,0) = to fix the coordinate freedom 
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C. Conditions on perturbations 

As in the case of self-similar solutions, we require (i) 
that the perturbations are analytic for all x G R, and (ii) 
that the perturbed space-times are regular at the center 
{Ap = at a; = — oo). 



D. The sonic point 

It can be seen that apart from the overall multiplica- 
tive factor, the perturbation solutions which satisfy the 
analyticity condition (i) at the sonic point are specified 
by one free parameter k, as we now explain. 

(1) First we note that the sonic point is a regular sin- 
gular point for the perturba tions . That is, the system 
of ode's for perturbations ( 5.13 ) is sing ular where the 
determinant of the matrix on the LHS of ( 5.13 ) vanishes: 
AxDx — BxCx = 0; this is identical wit h the sonic point 
condition for self-similar solutions, ( |4.5| ). 

(2) Next, in order to have a smooth s oluti on at the 
sonic point, the third and fourth row of ( 5.13 ) must be 
proportio nal t o each other, like in the case of self-similar 
solution, ( |4.6D . This yields an algebraic relation between 
^p(O), ^p(O), (I'p(O) and Vp(0) at the sonic point. 

(3) Third, we have an algebraic identity between 
Ap{x), Np{x), iDp{x) and Vp{x), given by ( ^.15| ). 

(4) The above two algebraic relations, together with 
our choice of gauge Np{0) = enables us to express, e.g. 
(Dp(0) and Vp{0) in terms of Ap(0); and then higher order 
expansion coefficients of perturbations are given in terms 
of ^p(O). Because the system is linear and homogeneous, 
overall normalization is irrelevant. We thus see that the 
solution which satisfies (i) is characterized by a single 
parameter k. 

This, together with the regularity condition (ii) at the 
center, in general allows only discrete values for k as now 
be explained. 



E. Asymptotic behaviour as a; — > ±cxd 

Searching numerically for eigenmodes of perturbations 
is in some sense easier than searching for self-similar so- 
lutions, because the ODE's are now linear and homoge- 
neous in unknown functions. In particular, asymptotic 
behaviour of solutions as a; ^ ±oo is easily obtained as 
follows. 



1. Asymptotic behaviour as x 



Substituting the asymptotic behaviour, (|4.12| ), of self- 
similar solutio ns as x — oo, the system of ODE's for 
perturbations ( 5.13| ) take on a simple form, given by 



Np 

LOp 



-1 
1 

2-3t^ 











-2/ 



Ap 



V„ 



(5.16) 



This homogeneous linear ODE with constant coefficients 
has four independent solutions, given by 



/ 



1 

-1 

2-37 
2(7-1) 





(5.17) 

The last one should be excluded in view of the identity 
( |5.15 ), and thus asymptotic behaviour of arbitrary solu- 
tions of ODE's is described by suitable linear combina- 
tions of the first three independent solutions. Now, one 
of them (the third one) blows up as a; — > — cxd. Therefore 
we see that any solution, which satisfies our boundary 
condition, should thus be obtained by choosing k so as 
to eliminate the unwanted expanding mode. 



2. Asymptotic behaviour as x —* oo 

Similarly, the asymptotic behaviour as a; — > cxd is ob- 
tained. It is now seen that there are three indep enden t 
solutions [after taking into account the identity (5.15)], 
whose a;-dependence are 



1, 



(5.18) 



as a; — + oo. It is thus seen that every solution stays 
bounded as long as Rck > 0, whereas for Rck < it 
blows up as a; — !■ oo. 



F. The gauge mode 

In searching for possible eigenmodes, one has to pay 
attention to a "gauge mode" which emerges from a coor- 
dinate transformation applied to the self-similar solution, 
as we now explain. 

Suppose one has a self-similar solution h(s,x) — 
Hss{x). Seen from anothe r co ordinate [related by the co- 
ordinate transformation ( |3.8D ], its perturbation is given 
by setting hp = in h of (^.12[): 



^var('5, x) 



{h = A,u;,V) 



(5.19) 



For general /i, this does not behave like an eigenmode. 
However, with the choice of /i(s) = e'"' {R arbitrary), h 
does behave as an eigenmode with eigenvalue k: 



3 (s, a;; k) 



A,Ld,V) 
N) 



(5.20) 
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The mode emerges via the coordinate transformation of 
the self-similar solution /igs, and thus can be considered 
as a result of pure gauge transformation, like a "zero 
mode" in translation invariant systems. Because k S C 
is arbitrary, this pure gauge mode forms a one-parameter 
family. 

It can also be seen that other (non-gauge) modes are 
transformed under the gauge transformation as follows: 
Eq.(5.19) for hvai{s,x) — hp(x)e'^'^ becomes 



{h^A,Oj,V) 
hp{x)e-^ + his)K,ix) + Ms) 
{h = N) 



As long as hp ^ 0, by taking /i(s) 



(5.21) 



this means 



transformed h is also an eigenmode, with the same k. 
Thus we have a one-parameter family of eigenmodes 
(with the same eigen value k) , which are mutually related 
by the gauge mode (5.20). 




FIG. 7. Profile of the eigenmode with the largest eigen- 
value. Curves represent Ap (solid line) , Np (dashed) , tjp (dot- 
ted), and Vp (dot-dashed). 



G. Numerical Results 



3. Other modes for ttie Evans-Coleman solution 



Based on the above observations, we searched for de- 
sired eigenmodes as follows. We first fix a value of k, 
and then starting from x = (the sonic point), integrate 
~oo. When there appear nonzero compo- 



( p^ to X 



nents of the expanding modes (for sufficiently large \x\), 
we judge that this k does not give a desired eigenmode 
and stopped. 

In this way, we searched for eigenmodes in the range of 
RcK > —1.5, |ImK| < 14. We found several eigenmodes 
which are explained below. 



1. The relevant mode for the Evans-Coleman solution 

The profile of the eigenmode with the largest Rbk 
is shown in Fig. |^. It has the eigenvalue n ~ 
2.81055255; which corresponds to the exponent value 
/3bh — .35580192, according to our scenario of sec- 
and which we believe to be exact to the last 



tion 
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To confirm our scenario, we performed a thorough 
search of other eigenmodes in the region < Re k < 
15, |ImK| < 14, and found none, except for the above 
mentioned relevant mode and the gauge mode. We also 
performed a less complete search in the region —1.5 < 
RcK < 0, |ImK| < 2. There is an eigenmode with 
RcK < —1.4, which is consistent with the results of the 
Lyapunov analysis in Sec. VI. 

It is not a priori obvious whether the eigenmodes form 
a complete set of basis functions. Due to the compli- 
cated structure of the equations of motion, we have not 
found a beautiful argument which can restrict possible 
eigenvalues (like that of |3^ and references therein). 

Our search of eigenmodes has a serious drawback that 
it is theoretically impossible to cover the whole values of 
K G C (unless, of course, one employs more sophisticated 
mathematical techniques). To fill this gap, and to further 
confirm our scenario, we performed a Lyapunov analysis, 
as explained in Sec. 0. 



digit. This is in good agreement with the value of g]. 



2. The gauge mode 



In our gauge, where Np{s, 0) = 0, we observe the gauge 
mode with R = —-^2^(0) ~ 0.35699, as explained in sec- 
tion ^[]. 



We take this opportunity to clarify the confusion expressed 
in "Note Added" of |ll|], which states that the gauge mode 



reported in (same as reported here) "seems to be erro- 
neous." The fact is that our report in jlo| (and of course 
here) is correct. The confusion seems to be due to different 
gauges used in the analysis. We here and in use the gauge 
Nyjs, 0) = at the sonic point, while Maison and we in 
Sec. use the gauge Np(s, — oo) = at the origin. The for- 
mer gauge gives the gauge mode at k ~ 0.35699, while the 
latter gives at k = 1; t he di fference in the values of k is well 
understood in view of (5.20). 
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H. Modes for other self-similar solutions 



A. Numerical Methods 



As has been stated in Sec. [V F , there exists a series of 



self-similar solutions (specified by the number of zeros of 
V) , in which the Evans-Coleman solution can be consid- 
ered as the first one (exactly one zero for V). We have 
searched relevant eigenmodes for the first several self- 
similar solutions and found more than one relevant modes 
(see Table |) . This implies that the other self-similar so- 
lutions are irrelevant for the generic critical behaviour. 
However, we here admit that our analysis is less com- 
plete for these higher self-similar solutions, because due 
to a difficulty in coding, (1) we have not done the analy- 
sis for (possible) self-similar solutions with even number 
of zeros of V, (2) there can be more relevant modes than 
is reported in the Table. 



self-similar solution 


K 


1 (Evans-Coleman) 


2.8105525488 


3 


8.456 




3.464 




1.665 




0.497 


5 


15.80 




7.13 




3.22 




1.51 




0.500 


7 


15.97 




6.92 




3.20 




1.50 



TABLE I. Relevant modes found for other self-similar so- 
lutions, which are labeled by the number of zeros of V. 



The Lyapunov analysis is a method of extracting eigen- 
values in descending order. It involves time integration 
the linearized EOM (^^) about the self-similar solution. 
It takes advantage of the fact that the eigenmodes with 
large Rck dominates at late times and that the volume, 
defined by an (arbitrarily chosen) inner product on F, 
of the parallelepiped spanned by the eigenmodes cor- 
responding Ki, K2,---, Kn approaches e'^i''2---K,i_ During 
the integration, vectors in F are orthonormalized by the 
Gram-Schmidt procedure. This is essential for avoid- 
ing numerical overflows and underflows. For details see 
Appendix |d|. We can also have information on the imag- 
inary part of eigenvalues by the period of oscillation in 
time s. 

We have performed the above-described calculation by 
numerically solving the PDE's for the perturbation, (^.4|). 
We e mphasize that solving the PDE's for perturbations, 
( ^.4|) , is much easier than solving the full EOM close to 
the critical point, because perturbations are well behaved 
in contrast with near-critical solutions whose derivatives 
are diverging in [t, r) coordinate. 

We used the first order Lax scheme, and the second 
order Lax-Wendroff scheme for time evolution. We show 
the results from the Lax-Wendroff scheme. To simplify 
coding, we employed a new gauge, iVvar(s, — oo) = 0. In 
this gauge we observe the gauge mode at k = 1. We have 
made appropriate coordinate transformations to stabilize 
the integration schemes. 

We imposed the free boundary condition at the sonic 
point ^ = 1. This can be done because no information 
can come in from outside the sonic point (as is easily 
seen by studying the characteristic curves), for spheri- 
cally symmetric metric perturbations contains no gravi- 
tational waves and are determined by the matter degrees 
of freedom. 



VI. LYAPUNOV ANALYSIS 

To further confirm the uniqueness of the relevant mode 
we performed a Lyapunov analysis around the critical so- 
lution. A Lyapunov analysis and the shooting method 
for an ordinary differential equation adopted in |l0|| and 
in the previous section are complementary to each other 
for the following reasons: (1) The former extracts eigen- 
modes in the descending order of its real part, whereas 
the latter affords information only of a finite region of 
complex K plane. (2) The Lyapunov analysis is useful 
even if the (discrete) eigenmodes do not form a complete 
set. (3) It is easier in the latter than in the former to 
numerically obtain accurate values of the eigenvalues k, 
hence the critical exponent /3bh- We here present our 
method of analysis and its results. 



B. Numerical results 

Fig. ^ shows the dependence of Rck on 7. The largest 
eigenvalue is real and agrees well with the first eigen- 
value (which is real) found by solving two-point bound- 
ary problem of an ODE in Sec. The second largest 
eigenvalue is that of the gauge mode, which should be ex- 
actly equal to one. For 7 < 1.2 the third and the fourth 
largest eigenvalues are real, and the fifth is imaginary. 
For 7 = 1.3, 1.4 the third is real and the fourth and the 
fifth is imaginary, which are complex conjugate, those 
for 7 > 1.4 are complex. We can consider that there is 
a cross-over of Lyapunov exponents corresponding to the 
real eigenvalue and the complex pair between 7 = 1.2 and 
7 — 1.3. For 7 > 1.5 the third and the fourth are a com- 
plex conjugate pair, so that there is another cross-over of 
Lyapunov exponents corresponding to the real eigenvalue 
and the complex pair between 7 — 1.4 and 7 = 1.5. For 
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to be discontinuous. (For the dust, 7 = 1, we expect a 
trivial critical behaviour with /3bh = 1-) This may be 
because the domain of attraction of the Evans-Coleman 
type self-similar solution vanishes as 7 — > 1+. 



FIG. 8. Dependence of Lyapunov exponents ReK on 7. The 
lines show the Lyapunov exponents. The top line, labeled 
"relevant" represents the relevant mode, and the second line, 
labeled "gauge," represents the gauge mode, whose Lyapunov 
exponent equals unity theoretically. For 7 < 1.3, the third 
and the fourth largest mode are real (labeled "real"). For 
7 = 1.3 and 7 = 1.4 the third is real and the fourth and 
the fifth are complex conjugate. For 7 > 1.4 the third and 
the fourth complex conjugate (labeled "cc"). We see that 
complex conjugate pair takes over the real eigenvalue twice in 
1 < 7 < 1.889. The graph shows there is a unique relevant 
mode for all values of 7 analyzed. 



large 7 it seems that many eigenmodes with close values 
of ReK are present. Our preliminary computation of ten 
modes for 7 = 1.889 shows that there at least nine modes 
in the range —0.6 < Re/t < —0.4, though it seems that 
not all of them are degenerate in ReK. 

The figure establishes the uniqueness of the relevant 
mode in collapse of perfect fluids with p — (7 — l)p, 
1 < 7 < 1.889. from this we can conclude that the critical 
behaviour must be observed for all of these models with 
1 < 7 < 1.889. 



VII. SUMMARY OF RESULTS 

We have performed the analysis presented in Sees. ^V| - 
VI for perfect fluid, with 1 < 7 < 2. As has been noted 
in EG a self-similar solution of Evans-Coleman type 



(i.e. with one zero of V) ceases to exist for 7 > 1.889. For 
Evans-Coleman type self-similar solutions, the Lyapunov 
analysis of Sec. ^ establishes that there is a unique rele- 
vant mode, and thus we can observe the critical behaviour 
(Fig. H). Then, the precise values of the relevant eigen- 
value K are obtained by the shooting method of Sec. 
and the critical exponent is given by /3bh = 

The result is shown in Table |l| This in particular 
shows that the relevant eigenmodes found by Maison | |Tl| ] 
are actually unique, and thus are responsible for the crit- 
ical behaviour. The value of the critical exponent /3bh 
depends strongly on 7. Moreover, the limit 7 — s- 1 seems 



7 


K 




/3bh 


1.00001 


9.4629170 




0.10567566 


1.0001 


9.45592488 




0.10575380 


1.001 


9.38660322 




0.10653481 


1.01 


8.74868715 




0.11430286 


1.03 


7.61774326 




0.13127247 


1.04 


7.16334221 




0.13959964 


1.05 


6.76491004 




0.14782163 


1.06 


6.41269915 




0.15594058 


1.08 


5.81789124 




0.17188358 


1.1 


5.33435815 




0.18746398 


1.12 


4.93282886 




0.20272343 


1.15 


4.44235059 




0.22510605 


1.18 


4.0484584 




0.2470076 


1.2 


3.82545008 




0.26140715 


1.22 


3.62729455 




0.27568756 


1.25 


3.36750228 




0.29695600 


1.28 


3.14337431 




0.31812947 


1.3 


3.00990875 




0.33223599 


1.32 


2.88714829 




0.34636253 


4/3 


2.81055255 




0.35580192 


1.36 


2.66838221 




0.37475891 


1.38 


2.57025726 




0.38906611 


1.4 


2.47850858 




0.40346844 


1.42 


2.39245265 




0.41798110 


1.44 


2.31150728 




0.43261815 


1.46 


2.23517329 




0.44739260 


1.48 


2.16301995 




0.46231659 


1.5 


2.09467339 




0.47740140 


1.52 


2.02980720 




0.49265763 


1.55 


1.93841621 




0.51588508 


1.58 


1.85338883 




0.53955219 


1.6 


1.79989076 




0.5555893 


1.62 


1.74873002 




0.5718436 


1.64 


1.69974510 




0.5883235 


1.66 


1.65278973 




0.6050376 


1.68 


1.60773076 




0.6219947 


1.7 


1.56444628 




0.6392038 


1.72 


1.52282404 




0.6566747 


1.74 


1.48276003 




0.6744180 


1.76 


1.44415717 




0.6924454 


1.78 


1.40692422 




0.7107703 


1.8 


1.37097467 




0.7294081 


1.88 


1.2383842 




0.8075039 


1.888 


1.2259859 




0.8156700 


1.889 


1.2244458 




0.8166960 




TABLE II. 


Values of n and /3bh 


= 1/k 


for 1 < 7 < 1.889. 


The last digit 


is rounded. 







VIII. CONCLUSIONS AND DISCUSSION 

The aim of this paper was twofold: First, in Sec. || we 
presented a mathematically rigorous framework, under 
which most of the critical behaviour observed in gravi- 
tational collapse are expected to be analyzed. We have 
in particular shown rigorous sufficient conditions under 
which (1) critical behaviour with a continuous self-similar 
solution is observed and the critical exponent is given ex- 
actly in terms of the un ique relevant Lyapunov exponent 
K by /3bh = (Sec. |II C| ), (2) critical behaviour with 
a discrete self-similar solution is observed and the crit- 
ical exponent is given by the same formula (Sec. |IID|), 
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and (3) different models exhibit the same critical be- 
haviour, i.e. univ ersality in the sense of statistical me- 
chanics (Sec. HE). The sufficient conditions for (1) and 
(2) are essentially the uniqueness of the relevant mode, 
together with some regularity on the dynamics. The con- 
dition for (3) is essentially that the difference of EOM is 
relatively small in the critical region. 

The key idea behind the analysis is the renormaliza- 
tion group. It naturally casts a system described by 
partial differential equations into a dynamical system on 
an infinite- dimensional phase space, where continuously 
and discretely self-similar solutions of the original equa- 
tions are described by fixed points and limit cycles. This 
opened the way to the abstract analysis on dynamical 
systems described in Sec. |l[ With the help of techniques 
developed in statistical mechanics to study critical be- 
haviour, we have been able to obtain detailed estimates 
on the time evolution of the deviation from the criti- 
cal solutions (or the stable manifold). We thus obtained 
general and precise results, including rigorous relations 
between the critical exponent and Lyapunov exponent. 

Another advantage of the renormalization group idea 
is that it enables us to treat systems with different EOM 
by naturally introd ucin g a space of dynamical systems, 
as described in Sec. II E . Each system is in general driven 
along a renormalization group flow in this space, and this 
enabled us to clarify terms of EOM which are "relevant" 
for the phenomena we are interested in. 

These point of view provide us a natural understanding 
on why critical behaviour are often observed in gravita- 
tional collapse. 

We then proceeded to confirm the pict ure for a con- 
crete model of perfect fluid in Sees. |lI-VII. We em- 
ployed essentially two methods to show that the main 
assumption, uniqueness of the relevant mode, is satisfied: 
(1) usual shooting method, and (2) Lyapunov analysis. 
Their results agree well, and strongly suggest that the 
relevant mode is unique. We o btained values of the crit- 
ical exponent /3bh in Sec. VII. 



needed only to confirm Assumptions Gl-3 in Sec. 11 C, 



It should be emphasized that the general framework 
presented in this paper is useful for both mathematical 
and numerical analysis on concrete models of critical be- 
haviour in gravitational collapse, because it tells that in 
which case one can rely on the intuitive picture of the 
phase space shown in the schematic diagrams Figs. |^, ^, 
^ and the argument presented in pO[ | which states that 
the linear perturbation is enough to understand all as- 
pects of the critical behaviour, apart from pathological 
behaviour mainly due to the infinite dimensionality of the 
phase space. For mathematical proof of the existence of 
the critical behaviour and estimation of the critical ex- 
ponent, the essential step is to find the norm || • ||, which 
defines the functional space or the phase space F, that 
satisfies the Assumptions in the Sec. II C . This should 



especially Assumption Gl. This is qualitative one rather 
than quantitative (because only the finiteness of the time 
before blowup is needed), and requires the time evolution 
of very small region of initial data, essentially a set of ini- 
tial data which is the fixed point plus the small relevant 
mode. The calculation of this e.g. in (s, x) coordinate 
should be much easier than adaptive mesh refinement cal- 
culation in (i, r) coordinates for wide region of the phase 
space and should allow a closer control of the numerical 
errors. This together with close analysis of the spectrum 
of the linear operator %j around the fixed point may lead 
to a computer-assisted proof of Theorem 

There remain several open questions to be answered. 
First, it sh ould be emphasized that our analysis in 
Sees. Ill -VII confirms only the local behaviour of the 



flow, around the specific self-similar solutions. Although 
self-similar solutions with more than one zeros of V seem 
to be irrelevant for generic critical behaviour (because 
they seem to have more than one relevant modes), it 
would be interesting to know the global behaviour of the 
flow, in particular, how these self-similar solutions are 
related. 

Second, now that various kinds of critical behaviour 
for different models have been observed, it would be in- 
teresting to ask what happens in a mixed system (e.g. 
perfect fluid -t- scalar flelds). We performed a prelim- 
inary analysis in this direction. That is, we consid- 
ered a "mixed" system of radiation fluid with a real 
scalar fleld, whose energy-momentum tensor is given by 

Tab = /5UaU6-|-p(UoWf,-|-5a&)+Va(/>Vb(/)-(l/2)5abV''0Vc'/), 

and studied the linear stability of the Evans-Coleman 
self-similar solution in this mixed system. (It is easily 
seen that eigenmodes decouple.) Lyapunov analysis in 
Sec. VI, applied to this mixed system, shows that the 



be highly nontrivial and depends on each model consid- 
ered. For numerical analysis the nonlinear simulation is 



Evans-Coleman self-similar solution has the unique rel- 
evant mode (given by the relevant mode of the radia- 
tion fluid), and all scalar eigenmodes are irrelevant. This 
means that the Evans-Coleman self-similar solution is 
stable under the scalar perturbation, and we will observe 
the same radiation fluid critical behaviour for the mixed 
system (at least when the scalar fleld is sufficiently weak). 
It would be interesting to investigate the behaviour of 
the system without the restriction that the scalar field 
is weak; in particular, where does the cross-over between 
two types (fluid, scalar) of critical behaviour occur? 

As is suggested from our general framework in Sec. ||, 
a "critical behaviour" such as observed in gravitational 
collapse is not in fact limited to self-gravitating systems. 
Similar, but slightly simpler, phenomena can be observed 
in much simpler systems, such as nonlinear heat equation, 
where we can rigorously carry out the analysis of Sec. ||. 
However if a solution of EOM blows up in such simple 
systems, it usually means that the equation is not physi- 
cally applicable in the blow-up region, and the blow-up is 
an artifact of bad approximation. General relativity pro- 
vides rare examples where the blow-up of solutions does 
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have a physical meani ng — for mation of a black hole 



Our analysis in Sees. [11- VII is not mathematically rig- 
orous. In view of the general framework given in Sec. O, 
and in view of the above remark that the critical be- 
haviour and blow-up of solutions are in fact physically 
meaningful, it is extremely desirable to develop a math- 
ematically rigorous analysis of critical behaviour for a 
physically interesting model such as radiation fluid or 
scalar fields. We are planning to come back to this prob- 
lem in the near future. 

Our intuition on gravitational collapse still seems to be 
heavily based on few exact solutions, especially the lim- 
iting case of pressureless matter. The critical behaviour 
may provide a different limiting case that the final mass 
is small compared to the initial mass for more realistic 
and wider range of matter contents. It will be of great 
help to settle the problems in gravitational collapse such 
as cosmic censorship conjecture. 



proves the continuity of Wr ■ Bounds on the norm follow 
immediately from the definition of ni. 

Step 2. Wr is " one- dimensional" , and any U, U' G Wr 
satisfy ||C/-C/'||, <C4|lC/-t/'||u. 

We first prove the claim for U, U' E Dr^n by induction 
in n, and then extend the result for general t/, U' G Wr- 
First, the claim is satisfied for n — Q, because Drfl = Dr 
is a straight line connecting C/q and Ui. Now suppose 
the claim holds for C/, U' G Dr.n and try to prove it for 
n„{U),'R.aiU') e Dr,n+i- [Any U G Dr,n+i should have 
at least one inverse image in Dr.n, so this is sufficient.] 
We c ompute as we did in deriving ( p. 8?: )~( 2.86D using 
( ^33|) , (write F = U' -U) 

n„{u') = n„{u + F) = na{u) + %,u'{F) 

+0{<oKi\\U -U*\\ \\F\\) + 0{Ki\\Ff). (A2) 

Now, by the inductive assumption, ||i^||s < H-F'llu = ll^ll- 
So we have 
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APPENDIX A: SKETCH OF THE PROOF OF 
LEMMA 2 

The proof consists of several steps, and is somewhat 
complicated. Because the Lemma itself is of auxiliary 
importance, and because the proof uses techniques which 
have been explained in detail in the proofs of other lem- 
mas and propositions already proved in Sec. 11 C , we here 
sketch only the essence of each step. 



Step 1. Wr is a continuous curve, connecting Uq and 

Un, + 1, With \\Uo - U*\\u = r, \\Un, + l - U*\\u > £2- 

We write Uq = U* + rF"'\ and C/„ = TlnaiU* +rF'''^). 
Because Uq satisfies the assumptions of Lemma |l|, 

• ||f/n — C/*||u keeps growing until it finally exceeds 
Si. This also shows ni{Uo; £2) < 00. 

• - U*\\s <\\Un- [/*||u for n < ni((7o; £2) + 1. 

Now Dr is a straight line connecting Uq and Ui ; and Wr is 
a union of images of Dr under TZna- for n < ni(t/o^£2) + l: 



Wr = ^nUDr.n, 



(Al) 



Because TZ^ is continuous and because ni is finite, TZna 
is also continuous (for n < ni). Thus Dr.n is a continu- 
ous set. Because endpoints of Dr^n are overlapping, this 



\n,iu')-n,iu)\u<e-^''\\F'\\ 

+Ki[6\\U-U*\\ + \\Fm \\F\\^ 

|7^,(c/')-7^,(c/)||„ 



(A3) 



> [e'^'^ - 6Ki\\U -U*\\- Ili^lU (A4) 



Taking \\U -U* 
leads to 



< £2 and < £2 into account, this 



\\n^iu')^n^{u)\\, 
\\n^{u')~nAu)\\u - e 



< 



\u-u% 



---7Ki£2 \\U-U'\\u 

7Xl£2 
' - £2 ■ 



(A5) 



It is now easily seen that this recursion preserves 

\\u-u'\\.. 



\U -U' 



(A6) 



and Cj < 1 follows from our choice of (52, £2, i-e. ( 2.65| ) 
and (|2.66| ). 

All this proves the claim for U, U' E Dr,n for some n. 
Because |jC/„ — U*\\ is strictly increasing, we can combine 
this bound for each Z)r,n, to conclude ||[/— /7'||s < C4IIC/— 
U'\\u for U,U' G Wr which are not necessarily in some 

Dr.n- 

Step 3. Contraction in -direction. We now consider 
the distance between a flow TZna{U) and Wr, and show 
that their distance in iJ'^-direction is decreasing. To be 
more precise, we introduce a decomposition (see Fig. S 
(c)): 



(A7) 



where this time C/„ G Wr{U*) and G £:"([/*). Because 
we have already shown that Wr{U*) is one-dimensional 
and extends from r to £2 in i5"-direction, we can always 
decompose this way, as long as r < ||[/("'^^ — J7*||u < £2- 
To derive a recursion for HG^H, we start as usual 
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= TZ„{Un)+%,u-.{Gl) + A (A8) 

where 

A = d{6K,\\U^ - U*\\ ■ \\Gl\\) + 0(A'i||G^||2). (A9) 

We have to define Un+i € Wr so that A is decomposed 

as 

A = Un+i-n,{Un) + F' {F'GE%U*)) (AlO) 

so that the entire expression can be written as C/„+i + 
Gl+i with 



Now ( |A10|) imphes 

l|A||u = ||^7„+l-7^,(^7„)||u 

and 

\\F^\\ = \\F%<\\Un+,-nAUn)\\s- 



lAI 



(All) 



(A12) 



(A13) 



However, we have shown in Step 2 that Wr is in some 
sense "parahel" to E'^iU*), i.e. 

\\u„+i ~ 7^,(^7„)||s < c4Un+i ~ n^{Un)\\u (au) 

(A15) 



Combining ( |A12| )-( |A14| ) imphes 

||Fi<(l + C4)|lA|| 
and we finally have 

n|G„|| + (1 + ^4) X 
0{6Ki\\Un-U*\\ ■ \\G. 



-0{K,\\Gir 

(A16) 



By ta king \\ Un ~ U*\\ and ||G„|| sufficiently small as in 
( 2.65 ) and ( 2.66 ), we can bound the right hand side by 
e-^"'"||G|| with some R" > 0. 

This in particular shows that intersections of Wr and 
the plane U* + E^{U*) + aF^°^ for various r's are Cauchy 
sequences, and thus the limit r \ exists: existence of 

Finally, repeating the above argument replacing Wr 
by ([/*), we get the contraction around W^^{U*) as 
well. □ 



APPENDIX B: SOLVING RECURSIVE 
INEQUALITIES 

We here provide a lemma, which effectively solves re- 
cursive inequa lit ies enc ountered in Sec. |ll|, such as (2.87), 
( p8|) , ( |2.174| ), ( 12.1751 ). 



Lemma 6 Let A, A, C , S and e be positive constants 
satisfying 1 < A^^ < A, G > 1, and 



e < 



A- 1 
~2C~' 



Let bn, fn be nonnegative sequences satisfying 
0<bo<e, /o = 0, 

bn+i iAb„±C [gMn V &„) + (/„ V 6„)'] 

fn+l <Ua+C [gnifn V &„) + (/„ V b„)^] , 

where gn (n>0) satisfies 

< 5n < S. 

Define n2 = maxjn € Z| 6^ < e for fc < n}. 

(i) n2 is finite and the following hold for n < • 



< /„ < b, 
K < e 



(a) Suppose in addition that gn satisfies 

gn < OO. 

n—0 

Then the following hold for n < n2 

< o«2 < 



3A- 1 
^ < feoA"^ < G'e, 

where C is some positive constant. 



(Bl) 



(B2) 
(B3) 
(B4) 



(B5) 

(B6) 
(B7) 

(B8) 

(B9) 
(BIO) 



The proof is an elementary exercise in calculus, which 
is done in several steps in a bootstrapping manner. We 
only consider n < n2, which in particular means < 
b„ < e. 

(i) Step L Rough estimates and f„ < b„. We first 
prove by induction /„ < 6„ (i.e. (B6)), and 



A + 1 , , 3A - 1 

r 6„ < bn+l < b„ 



(Bll) 



These are satisfied for n = 0, because fo — 0, and be- 
cause of our choice of e. Now suppose they hold for n. 
Our inductive assumption in particular includes /n < &n, 
which reduces the recursion (B3) to 



bn+l ^ Abn ± Cbnign + bn) 

^ 5„[A±G(<5 + e)] 



(B12) 
(B13) 



where we used uniform bounds on gn and 6„. Eq.(Bll) 
immediately follows from our choice of S and e. 



On the other hand, /„ < 6„ reduces the recursion (B4) 



to 
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fn+1 < A/„ + Cbn{gn + b„ 

< 6„ [A + CiS + e)] 

< [A + CiS + e)] 



A + 1 



(B14) 



where in the last step we used ( Bll ) just proved. The 
f act or before 6„+i is easily seen to be less than one, and 
(B6) is proved. 



Step 2. Rough estimate on bn- First note that (Bll) 
for n = n2 and the definition of n2 implies 



e < b. 



< 



3A- 1, 



(B15) 



i.e. (B£). Also solving (Bll) backwards from n2 to n, 
we have 



bn2 



< 



A + 1 



(B16) 



which, in view of < e, is nothing but (B7). 

Step 3. Refined estimate on 6„. For later use, we note 
the following refined estimate, which is valid without the 
additional assumption ( p38| ). First note that g„ < (5, < 
e, and /„ < &„ gu arant ee K-C{gk + bk) > (A + l)/2 > 1. 
So the recursion (B12) implies {n < 712) 



ri-l 



n [A±C{gk + b,)] 



k=0 



< 



fe=0 



1 ± 



C{gk + bk) 



A 



§ A" exp 



± 



(21n2)C 
A 



k=0 



(B17) 



In the second step above we used | ln(l + x)| < (2 ln2)|a;|, 
which is vahd for < 1/2. [Note that C{gk + bk)/A < 
(A-l)/(2A)<l/2.] 

(ii) Given the above estimat e (iBlTl ), the rest is easy. 
Making use of the assumption ( |B8D on g„, and the esti- 
mate (B7) on bk proven so fa r, it can be easily seen that 
the sum in the exponent of (B17) is bounded by a con- 
stant depending only on C, A, A and XliLo 5^' uniformly 
i n n a nd n2 as long as n < n2. Taking (B£) into account, 
( ^17| ) with n^n2 gives ( pl0| ). □ 



APPENDIX C: SPHERICALLY SYMMETRIC, 
SELF-SIMILAR SPACE-TIME 

Defriese and Carter (see Q|) showed that the line el- 
ement of any self-similar, spherically symmetric space- 
time is written as 

ds^ = e-'^'^{-Gi{y)dTi^ + G2{y)dy'^ + e^y{de^ + smed(t)'^)). 

(CI) 



Let t 



-W{y) and 
Y'{y) = / 



e'^ where 



dy- 



J constant 



2G2 



One has 

ds^ = -a^dt^ + a^dr^ + r^{d6'^ + i 



(C2) 



(C3) 



where a and a are functions of y. Since r/ (— t) is a func- 
tion of y only, a and a are functions of r/{—t). 



APPENDIX D: THE LYAPUNOV ANALYSIS 

The Lyapunov method has been extensively used in 
nonlinear analysis |^5| . To make our presentation self- 
contained and transparent, we present a concise mathe- 
matical description of the method, restricted to the situ- 
ation we are interested in. 

Let us consider the linearization of a flow on a real 
Hilbert space (or a real Hilbert manifold) around an or- 
bit. It is determined by an equation 



ds 



{s)=Ais)f{s), 



(Dl) 



where each /(s) is an element of real Hilbert space V, a 
complete vector space with inner product, and A{s) is a 
real function determined by the (background) orbit which 
we are considering. Let (•,•) denote the inner product 
and II • II = denote the norm defined by the inner 

product. We have A{s) = A in ( pl[ ) when the orbit is a 
fixed point, and A{s -|- A) = A{s) when it is an periodic 
orbit with periodicity A. We concentrate on these cases 
below. 

Let us define time evolution operator : ^ by 



/(O) = F, fis) = T,F, 



(D2) 



where / is a solution of (Dl). We wish to find eigenvalues 
K and eigenvectors E"^ of Ts satisfying 



(D3) 



in particular those with large Rck. Here k is a complex 
number and E"^ is a complex vector with KeE"^ ,lmE'^ g 
V. In (D3) and below, s (and a) can be any number in 
the case of a fixed point, whereas it is restricted to be 
integer multiple of A in the case of a periodic orbit. 

We assume that the eigenvalue problem ( p3| ) has 
a complete set of discrete eigenvectors, denoted by 
S'^ = {-Ej'^}i=i,2,..., with the corresponding eigenvalues 
{Ki}i=i.2...., RcKi > ReK2 > ReK3 > .... Each Rck^ is 
called the ith Lyapunov exponent. (Here we implicitly 
assume that we can always find Ki with the largest real 
part after we have defined k^-i.) We observe that 

if Ki is not real there is an integer i* ^ i with Ki* = m* 
and Ef, = Ef*, where asterisks (except for those on i or 
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j) denote complex conjugate in this appendix. We define 
£ = {-Ei}i=i,2,... by Ei = E'[ if is real, and Ei = ReE'^, 
Ei* = Imi?j'= if Ki is not real, i < i* and E^, — E". Any 
real vector F can be expanded by £ or £'^ as 



(D4) 



where = /f * = (1/2) (/' + ^/^^'). Without loss 
of generality, we assume all eigenvectors are normalized 
so that ||Re£;,'=f + ||Im£;,'=||2 = 1, although E,'s are not 
necessarily orthogonal with respect to the inner product 

(•,•)• 

Let J- — {i^i}i<i<„ denote an n-frame of (linearly in- 
dependent, real) vectors and let spanJF be the subspace 
of V spanned by J^. Below we often omit 1 < i < n in 
T — {^i}i<i<n if no confusion occurs. 

We define several operations on frames. First given an 
operator X : V t—^ V we define X, an operator on frames 
induced from X, by XT = {X{Fi)}i. We in the following 
consider the cases where X is a linear operator (such as 
To-), or the normalization operator A'': NF = 
Second, we define the orthogonalization operator OJ- t-^ 
T' — {i^/}i<i<n, where T' is defined by 



||P'||2 
.7 = 1 " 



2<i<n. 



(D5) 



We adopt the convention that XTi = {X!F)i denote the 
ith vector of n-frame XT , when X is one of the opera- 
tions on a frame defined above. Distinct from X induced 
from X : V V , XTi<i<n in general depends on Fj 
with i ^ i- Note that spanOJF — spanjT. 



Proposition. For almost every T , i.e., except for "mea 
sure zero" cases whose precise meaning is given by (D16) 
we have 



r 1 , \P%naJ'^\ 

iieKi = lim In — 

m^oo ma Fj 



^ m 

lim y2Xi{l,(j), 



va^ca ma 



1=1 



for 1 < i < n, where 

K{t,a) = \n\\0%{NO%)'-^T\ 



(D6) 
(D7) 

(D8) 



Remark. The first expression (D6) of the Proposition 
is mathematically straightforward, but is not suitable for 
numerical calculations, due to serious overflow and un- 
derflow problems in large time integration. To overcome 
this difficulty, we employ the second expression of the 
Proposition, which is performed in the following manner 
in practice. 



(0) Prepare T and let = 0. 

(1) Evolve T in time by T^, and define T' = T^rT. 

(2) Find T" = OT„T and P" = MOT^T by the 
Gram-Schmidt procedure: 

F1 = F[^{F';',F-)F';\ F^' ^F':ilF':i, 

i > 2. (D9) 

(3) Add A, = In II F/' 1 1 to Ai. 

(4) Define new T = T'" and go back to (1). 

XJma, where m is the number of iteration, gives the 
Lyapunov exponent Rck^ according to the Proposition, 
because NO is equivalent to the Gram-Schmidt proce- 
dure. 

Before proceeding to the proof of the Proposition we 
list several properties of operations on frames. 

Lemma. Let X ,y be operations on frames. 

(i) If X is induced from a linear operator X : V ^ 
V , then as operators on frames, OXO = OX . In 
particular. 



O^ = O, o%o = o%. 



(DIO) 



(ii) If X,y are induced from operators X,Y on V, we 
have Xyj- = {X {Y {Fi))} In particular, 

N'^^N, N%N = N%. (Dll) 

%NT = {\\F.,\\-^Ts{Fi)}i<,<n. (D12) 



(Hi) For Ci ^ (I < i < n) we have 

0{{c^F,}l<^<n) = {c^OT^}l<^<n (D13) 

and 

NON^NO. (D14) 

Proof. (i) By the definition of O, any frame T is 
written as = [OF^ + Y.'~X ahOFj] with suit- 

able scalars a^iS. Since X is linear we have XT = 
I XOTi + y\l~\ a^iXOTn } . However, for scalars 

L J l<i<n 

O { \F, + y^a^,FA I =0({Fji<,<„), (D15) 

J = l 
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which can be shown by induction in i. So we have 
OXT = OXOT. 

(ii) is just a definition. 

(iii) can be shown by induction in n, noting that OTi 
depends only on Fj (j <i). □ 

Proof of the Proposition. We p rove the proposition 
in two steps. First, we prove (Dq), an d th en prove the 
equivalence of expressions in ( |D6D and (D7). 

We assume T to be generic in the sense that 

span{i^j}j<i n (span{£'j}j>j\span{£'j}j>i+i ^ {0}), 

l<i<n. (D16) 

The set of all nongeneric J^'s is of measure zero in the 
space of n-frames, because each nongeneric T contains 
at least one such i that Fi lies in a subspace of V of 
codimension one, spanji^jlj^j. from ( D16| ) we can define 
= {Fi\i such that (for \ < i <n) 

Fi e span{£'j}j>i\span{£'j}j>j+i, 
i^,-F,Gspan{F,}j<,_i. (D17) 

Note that we have 

OT = OT. (D18) 
Let us expand T in terms of E as 



Oi ^ 5 ^ O2 



(D24) 



for sufficientl y lar ge s, where Ci, C 2 ar e positive con- 
stants. Eqs. (D26|) and (|D24D prove (ra). 



Now we tu rn to the proof of equivalence of (D6) and 
Using (|D11|) and (|D14|) we have 



A/'OTAA = NOAfTN = NONT = NOT 
and thus 

(AAOT)^ = N[0TY. 
So by the definition of cr), we have 

A,(/, a) - In \\Or„N{OT,)'-^T,\\ 



= In- 



(Or,)'"l.F,;|l' 



(D25) 



where in the last step we used (D12) and (D13). We thus 
have 



^ Ai(Z,CT) = In 



(D26) 



and the proposition is proved. 



□ 



wher e /^-, 's are defined similarly as /"s in (D4). Condi- 
tion ( D16 ) implies that /% 7^ for all i < n, and that no 
/■'i's appear for j < i. Time evolution of Fi is given by 



00 

r,F, = ReJ^e'^'f^E^ = e^°'''^ + 0{e-'' 



where 

Ei{s) = Re 



E 

j: j>'i,RcKj— Rc/t; 



(D20) 

g»(Im«,).^j^^c^ (D21) 



and (5 is some positive constant. Since 7^ 0, ||i?i(s)| 
is bounded and nonzero uniformly in s. Operating O on 
both sides of (D21), we thus have by (D13), 



(oi,{s) + Oie-'')) 



(D22) 



for sufficiently large s. The left hand side can be replaced 
by OT sTi leading to 



OTsT, = e 



RcK.j S 



(0E,{s) + 0{e-'') 



(D23) 



because we have from ( DIP ) and ( D18 ) that OT sT = 
OTsOT = OT.OT = OT^T for any s. Thus we have 
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